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Abstract 

In the present article we study the form factors of quantum integrable lattice 
models solvable by the separation of variables (SoV) method. It was recently 
shown that these models admit universal determinant representations for the 
scalar products of the so-called separate states (a class which includes in partic¬ 
ular all the eigenstates of the transfer matrix). These results permit to obtain 
simple expressions for the matrix elements of local operators (form factors). 
However, these representations have been obtained up to now only for the com¬ 
pletely inhomogeneous versions of the lattice models considered. In this article 
we give a simple algebraic procedure to rewrite the scalar products (and hence 
the form factors) for the SoV related models as Izergin or Slavnov type de¬ 
terminants. This new form leads to simple expressions for the form factors in 
the homogeneous and thermodynamic limits. To make the presentation of our 
method clear, we have chosen to explain it first for the simple case of the XXX 
Heisenberg chain with anti-periodic boundary conditions. We would neverthe¬ 
less like to stress that the approach presented in this article applies as well to a 
wide range of models solved in the SoV framework. 
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1 Introduction 


Quantum integrable systems are ubiquitous in modern theoretical physics appearing 
both in statistical mechanics and field theory with applications ranging from con¬ 
densed matter to string theory [MQ]. They provide unique possibility to obtain 
non-perturbative and exact results for strongly correlated systems that cannot be 
obtained by other methods. Besides the computation of spectrum, scattering ma¬ 
trices and partition functions, one of the main challenges in this domain concerns 
the exact computation of the form factors and correlation functions that connect to 
measurable physical quantities in these systems. 

The quantum inverse scattering method [nHn together with its associated 
Yang-Baxter OESHIZ] and quantum group structures [I8ll22] provide a powerful 
framework to tackle such problems. A central object in this approach is played by 
the so-called quantum monodromy matrix T(A) depending on a continuous complex 
parameter A whose matrix elements are operators acting on the quantum space of 
states of the systems of interest V and satisfy quadratic commutation relations gov¬ 
erned by an R matrix solving the Yang-Baxter cubic equation. The main point of the 
method is the existence of an abelian sub-algebra of operators acting on V, including 
the Hamiltonian of the system, generated by the transfer matrix T(A) constructed 
algebraically from T{X) and leading, through an expansion in A, to a complete (in the 
sense of characterization of eigenstates) set of conserved operators and (dynamical) 
symmetries responsible for the integrability of the model at hand. In this framework, 
the solution of the original spectral problem for the Hamiltonian is embedded in the 
A-independent resolution of the spectral problem for the transfer matrix T(A) as a 
linear operator on V. At this point, within these algebraic settings, several methods 
can be used to solve this spectral problem. The first methods historically, Bethe 
ansatz [T] and algebraic Bethe ansatz mHH], have been extensively used with great 
success in many paradigmatic integrable systems like Heisenberg spin chains or sev¬ 
eral lattice discretization of integrable field theories like sine-Gordon. It appeared 
however that for systems lacking an obvious reference state, namely a simple eigen¬ 
state of the transfer matrix T(A) from which the complete space of states V can be 
generated by successive actions of matrix elements of the monodromy matrix T{X) 
on it, a different, somehow more generic, approach should be designed. The quantum 
separation of variable (SoV) method initiated by Sklyanin [231126] provides such a 
powerful tool for computing the spectrum and the eigenstates of quantum integrable 
systems, especially for systems lacking such a reference state. Sklyanin pioneering 
works on SoV have been presented for fundamental examples of integrable quantum 
lattice models, e.g. the spin 1/2 XXX chains and the Toda chains, after them sev¬ 
eral contributions have brought to further develop of the SoV method and nowadays 
several classes of integrable quantum models are proven to admit a description in 
the framework of the Sklyanin’s SoV method [27H5I]. One very important feature of 
the SoV method is that it provides not only the equations determining the spectrum 
of the transfer matrix but also the proof of its completeness and the construction 
of the corresponding eigenstates; this is in contrast to the Bethe ansatz or algebraic 
Bethe ansatz approach where the proof of completeness is in general a non trivial 
task [52l[53]. 

For all these models, one is of course interested to go beyond the knowledge of 
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the spectrum properties so as to reach their dynamical behaviour through the com¬ 
putation of their form factors (matrix elements of local operators in the eigenstates 
basis of the transfer matrix T(A)) and correlation functions. Several progresses in 
this direction have been achieved in the recent years, in particular for solvable lattice 
models |54H64j . As a matter of fact, the computation of form factors and correlation 
functions needs an explicit representation of the local operators in the eigenstate 
basis of the transfer matrix T(A) together with manageable formulas for the com¬ 
putation of the resulting scalar products of states. While the first problem can be 
solved within the quantum inverse scattering method by essentially expressing local 
operators in terms of the quantum monodromy matrix entries, hence solving effec¬ 
tively the so-called quantum inverse scattering problem |571I65| . the answer to the 
second question was up to now strongly dependent on the method used to describe 
the spectrum and eigenstates of the system. 

In the context of algebraic Bethe ansatz the computation of form factors and 
correlation functions [57ll59ll65l[67] was performed in settings sufficiently explicit to 
compute their critical behaviour [6814721 and make explicit contact (at operator level) 
with conformal field theories m- All this program was finally pushed forward to 
time and temperature dependent correlation functions [Z1H79], not mentioning clear 
contact with structure factors accessible in particular through neutron scattering 
experiments on magnetic crystals [8Ull81j . An essential technical feature in all these 
results was the appearance of rather sophisticated determinant formulas, in partic¬ 
ular for representations of certain partition functions [52], scalar products of states 
and form factors [571183] , in a form suitable for their analysis in the thermodynamical 
limit one is usually interested in [684F7TlF77ll841[85] . 

In the SoV approach, this programme is up to now slightly less developed due 
to peculiar technical features. To explain this point in more detail, let us recall 
that in the SoV approach one needs first to identify a diagonalizable operator, say 
S{X), computed from the monodromy matrix T{X) entries, having simple spectrum, 
in such a way that the multidimensional spectral problem for the transfer matrix 
T(A) separate in the 5(A) eigenstate basis into multiple one dimensional spectral 
problems, hence leading to its solution. In the example of one dimensional lattice 
models, the space of states V is realized as a tensor product of local quantum space 
of states Vn at sites n of the lattice, n = 1,... iV, V = ®n=i^n- In this case, in 
order to find a suitable operator 5(A) with simple spectrum, one is often led to 
introduce inhomogeneity parameters attached to each site n of the lattice and 
to let them be in generic positions in the complex plane. This has no effect on the 
integrability properties of the model and this is somehow a standard method used 
in integrable systems; namely to solve a given model by embedding it in a larger 
class of models depending on continuous parameters, here the and use these new 
variables to obtain simpler resolution process. All computations are then performed 
for this inhomogeneous integrable model, the homogeneous limit being taken only 
at the very end on the quantities one is interested in. It has been first shown in [86| 
that within such an SoV framework, the scalar product of separate states and even 
the form factors of local operators admit simple determinant representations written 
in particular in terms of Baxter Q functions. The key role in these representations 
is played by a Vandermonde determinant and its various (straightforward) dressing, 
making the computation of scalar products in the SoV framework much simpler and 
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transparent compared to the one obtained in the algebraic Bethe ansatz settings. 
Such representations have quite universal character and properties as it has been 
shown in many following examples |44ll46lH8ll87ll88j . Moreover it can be anticipated 
that such features will generalize to even more complicated systems associated to 
higher rank quantum groups, making the SoV approach even more appealing. The 
power and applicability of all this method however requires that the determinant 
representations for the quantities like the form factors or correlation functions are 
simple or explicit enough in terms of the inhomogeneity parameters such that 
the required homogeneous limits are non ambiguous. In particular one would like to 
avoid implicit zero over zero expressions in the homogeneous limit that could just 
spoil the simplicity and benefits obtained within the SoV method. However, it has 
been observed already in [86ll87j . that the homogeneous limit of the determinants 
for scalar products or form factors of local operators obtained by suitable dressings 
of the Vandermonde determinant do not have obvious homogeneous limit although 
there is no doubt that such a limit should, by construction, be well defined. 

It is the purpose of the present article to show that the rather universal determi¬ 
nant representations of scalar products of separate states obtained in the framework 
of SoV can be generically recast in a different form allowing for an obvious ho¬ 
mogeneous limit, hence opening the way to use extensively this approach to tackle 
dynamical properties of quantum integrable systems. More precisely, it will be shown 
that the dressed Vandermonde determinants obtained in the SoV approach for scalar 
products are equal to certain Izergin-type determinants, making the above discussed 
homogeneous limits trivial. These representations are shown in their turn to be 
equal to generalized Slavnov-type determinants if one of the two separate states is 
defined by a solution to equations of Bethe ansatz type. All these equalities between 
determinants are obtained through purely algebraic identities, hence making their 
appearance quite universal. For simplicity and also for making the correspondence 
between SoV results and algebraic Bethe ansatz ones explicit, we made the choice 
to explain the essential features of our results in a very elementary example : the 
XXX anti-periodic chain, in this case these formulae are reminiscent of those ob¬ 
tained in [5211221 for the XXX periodic chain. Recently several authors studying 
in particular the properties of the Slavnov formula for the XXX model discovered a 
relation of this type of determinant with Vandermonde type representations mm- 
There are also some more complicated relations observed for the XXZ chain |93] . 
We would like to stress however that the method presented in this paper seems quite 
universal leading to similar formulas for more complicated cases, like for XXZ or 
XYZ Heisenberg chains, even in the presence of generic integrable boundaries, that 
will be given in a separate publication. 

The article is organized as follows. In section 2 we recall the basics of the SoV 
method for the anti-periodic XXX model. In section 3 we describe the determinant 
formulas for scalar products of separate states that include all eigenstates of the 
transfer matrix. We then derive the main identities and equivalent representations 
of these scalar products in terms of Izergin-type and generalized Slavnov-type de¬ 
terminants when one of the state is an eigenstate of the transfer matrix. Similar 
formulas are then obtained for form factors of local operators in section 4. In section 
5 we show the equivalence of this SoV description with the one obtained from alge¬ 
braic Bethe ansatz, making use of the fact that by an explicit and simple change of 
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basis the anti-periodic XXX model can be recast in a form allowing its resolution 
by algebraic Bethe ansatz. 


2 The antiperiodic XXX Heisenberg chain in the SoV 
framework 

In this section we introduce the XXX Heisenberg chain with antiperiodic boundary 
conditions and recall the solution of the inhomogeneous version of this model by 
means of Sklyanin’s Separation of Variables approach |23ft26] . Within the frame¬ 
work of this approach, the eigenvalues and eigenstates of the transfer matrix are 
completely characterized by the solutions of a system of discrete equations involving 
the inhomogeneity parameters of the model, which can be understood as discrete 
versions of Baxter’s famous T-Q equation [3]. We show here that this SoV discrete 
characterization of the transfer matrix spectrum and eigenstates can be reformulated 
in terms of polynomial solutions of the continuous (i.e., functional) T-Q equation. 
This enables us to obtain a complete description of the spectrum in terms of the 
solutions of a system of Bethe-type equations, and to rewrite the corresponding 
eigenvectors in an ABA-type form, which we expect to be more convenient for the 
consideration of both the homogeneous and thermodynamic limits of the model. 


2.1 The antiperiodic XXX Heisenberg chain 

In this paper we consider the XXX Heisenberg chain of spin 1/2, 


N 

+ «+l - 1 ) > 

n=l 

with antiperiodic boundary conditions. 


0-^+1 = a = x,y,z. 


( 2 . 1 ) 


( 2 . 2 ) 


The Hamiltonian (12. ip acts on a 2'^-dimensional qnantum space V = with 

Vn — C^. Here and in the following, (t“, a = x,y,z, stand for the Pauli matrices at 
site n i.e., acting on the local quantum spin space Vn- 

The iZ-matrix of the model corresponds to the rational solution of the Yang- 
Baxter equation, 


R{X) 


/X+y 0 0 0 \ 

0 X r] 0 

0 T] X 0 

Vo 00 X+y / 


(2.3) 


where A is the so-called spectral parameter and y is an arbitrary non-zero complex 
parameter. The monodromy matrix of the XXX spin-1/2 chain is defined as the 
following ordered product of iZ-matrices, 


To{X) = iZo7v(A — Cn) ■ ■ ■ Roi{X — .^i) 


(AiX) B{X)\ 
\C{X) D{X)) 
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(2.4) 




This is an operator which acts on the tensor product of a two-dimensional auxiliary 
space Vb — with the 2^-dimensional quantum space V of the model. We have 
used here the standard notation which may label by indices (at least when it is not 
clear by the context) the space(s) on which the corresponding operator acts in a non¬ 
trivial way. Note that we have introduced in (12.41) a set of inhomogeneity parameters 
Cj) 1 ^ j ^ -^5 so that To (A) (12.4p corresponds in fact to the monodromy matrix 
of an inhomogeneous generalization of the XXX spin chain. These inhomogeneity 
parameters are crucial for the SoV study of the model. 

The transfer matrix of the model with antiperiodic boundary conditions can be 
defined as follows 

r(A) = tro (aS To(A}) = B{\) + C(A). (2.5) 

This equation defines a one-parameter family of commuting operators. It is impor¬ 
tant to underline that the transfer matrix is a polynomial in A of degree N — 1. 
It is also easy to observ^ that the transfer matrix satisfies the following simple 
symmetries, 

[5^r(A)] = o, [r^r(A)] = o, ( 2 . 6 ) 

where 

S- = al, r* = < = (-i)" exp (^sA . (2.7) 

In the homogeneous limit (^n 0), the Hamiltonian (|2.ip of the XXX spin quantum 
chain with antiperiodic boundary conditions is obtained as a logarithmic derivative 
of the antiperiodic transfer matrix (|2.5I) : 

(2.8) 

A=0 

Finally, let us recall the expression for the quantum determinant of the monodromy 
matrix, which is a central element of the Yang-Baxter algebra, 

det To(A) = B(\)C{\ - r?) - A{\)D(\ - t]) = -o(A) d{\ - rj), (2.9) 
<? 

with 

N N 

a{X) = ll{X-Cn + v), d{X) = lliX-Cn)- ( 2 . 10 ) 

n=l n=l 




2.2 Separation of variables for the XXX spin chain 

The functional version of the separation of variables solution of the spin-1/2 XXX 
spin chain leading to the diagonalization of the transfer matrix (|2.5p comes back 
to the early works of Sklyanin |251I26] . Following [46] . here we briefly recall the 
construction of the SoV basis of the space of states as well as the characterization 
of the transfer matrix spectrum and eigenstates which follow from this approach. 

Starting from this point we will always suppose that the inhomogeneity param¬ 
eters Ca, 1 < o < X, are generic, or at least that they satisfy the condition 

for /iG{0, 1 }, ya^b, (2.11) 

^See the proof of Theorem 14.21 for an explicit derivation of these symmetries. 
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which ensures the validity of the SoV approach (see [25l|26]). 

The SoV diagonalization of the transfer matrix (j2.5p relies on the construction of 
a basis of the space of states, that we shall call SoV basis, which for the antiperiodic 
case can be chosen as the one diagonalizing the action of the operator D{X) of the 
monodromy matrix. In this basis, the action of the operators B{X) and C(A), and 
hence of the antiperiodic transfer matrix (I2.5p . happens to be qnasi-local. In the 
case of the XXX spin-1/2 chain, the construction of such a basis can be explicitly 
realized in the quantum space V (respectively its dual space) by a multiple action of 
operators B{^j) on the reference state (respectively of operators C{^j) on the dual 
reference state). In fact, as noticed in this SoV basis is nothing else that the 
so-called F-basis introdnced in [95] . 

Concretely, let us define, for each V-tuple h = (hi,..., h^), a state | h) in the 
quantum space V and a state (h | in the dual quantum space as 


|h) 

(h| 


1 

1 




( 2 . 12 ) 

(2.13) 


where | 0) and (0 | stand respectively for the ferromagnetic reference state with all 
the spins up (right reference state) and for its dual state (left reference state): 


Af / 1 \ N 

|0) = ® , (0| = ® (1,0)[„]. (2.14) 

V ^ / [ n ] ’^=1 

For convenience, the normalization factor in (|2.12D and (|2.13D is chosen to be a 
Vandermonde determinant of the inhomogeneity parameters of the model: 


nm)= n (2.15) 

l<b<a<A^ 


The condition (j2.1ip on the inhomogeneity parameters of the model ensures that 
the operator D{X) is diagonalizable and has simple spectrum. It can easily be shown 
that the 2^ states | h) constitute a complete set of eigenstates for this operator: 


F(A)|h)=dh(A)|h), 


(2.16) 


and similarly, 


(h|F(A) = dh(A)(h|, 


(2.17) 


where 


N 

dhW = + Kr]). (2.18) 

n=l 


Moreover, the basis of the quantum space built from the vectors (|2.12|] and the 
basis of the dual quantum space built from the vectors (I2.13p are orthogonal. More 
precisely, with the chosen normalization, the scalar product of a state | h) of the 
form (|2.12p and of a state (k | of the form (|2.13p is given by 


_ (^k,h 

“ vm)v{{^-hr,}y 


(2.19) 
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where stands for the Vandermonde determinant (|2.15p associated with the 

set of inhomogeneity parameters ^j, 1 < j < N, whereas V ({^ — hr]}) stands for the 
Vandermonde determinant associated with the set of shifted inhomogeneity parame¬ 
ters — hjT], 1 < j < iV. Hence we have the following decomposition of the identity 
I on the quantum space V: 


^ = Vm) E V{{^-hr]})\h){h\. 
he{o,i}^ 


( 2 . 20 ) 


From the SoV basis (j2.12p or (j2.13p we define as in [IB] what we call separate 
states. These are states which admit a particular factorized form when expressed 
in these basis, so that the computation of their scalar product will follow quite 
straightforwardly from (I2.19p . More precisely, these are states which can be written 
in the form 


N 


^“ 1 = Yl(^i^a-har])V{{(-hr]}){h\, 

hg{0,l}^ a=l 


N 


1/3) = Yh^i^a-haT]) v{{^ + hr]})\h), 

hg{0,l}^ a=l 


( 2 . 21 ) 


( 2 . 22 ) 


for any function a or /3. It is important to mention here that such states are defined, 
up to a global normalization, only by the N ratios a{^j—r])/a{^j) (or 
1 < j < V. This means in particular that many different functions a may lead to 
the same separate state. 

For the effective computation of the scalar product between states of the form 
(|2.21l) and states of the form (I2.22p . it is convenient to rewrite (I2.22p with the same 
Vandermonde determinant as in (I2.2ip by means of the property 


N 


N 


(-vn n 


n=l 


Cn Cm T f] 


Cn Cm f] ^ 


Vi{^-hv}) = Vi{C + hr]}), 


(2.23) 


which can be proven by direct computation. Then the right separate state (j2.22l) 
can alternatively be written as 


N 


1/3) = 

hg{0,l}^ a=l 


(2.24) 


where the function /3 is such that 


^(Cn) = /3(Cn), mn-r]) = - 


a(Cr) 


d{in - v) 


/3(Cn - r]). 


(2.25) 


The eigenstates of the antiperiodic transfer matrix (|2.5I) happen to be particular 
cases of separate states, associated with some function (5(A) with ratios Q{Cj — 
'n)/Q{^j)j 1 < J < -/V) fixed by the corresponding eigenvalue of the transfer matrix. 
More precisely, we can formulate the following theorems, which were proved in 
and in |l6| (for a more general XXZ case). 
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Theorem 2.1. Let the inhomogeneity parameters satisfy the eondition (12.lip . Then 
the antiperiodic transfer matrix T (A) (12.51) has a simple spectrum, and a funetion 
r(A) is an eigenvalue ofT{X) if and only if 

1. it is a polynomial of degree N — 1, 

2. it satisfies the following set of quadratic discrete equations: 

r{f,n)T{Cn-v) + a{f,n)d{Cn-v) =0, Vn G {1,..., iV}. (2.26) 

Theorem 2.2. The left eigenstate of the transfer matrix T{X) corresponding to the 
eigenvalue t(A) is a left separate state {Q\ of the form (|2.2ip for a function Q{X) 
satisfying the following set of conditions for all n £ {1,..., N} : 

r(Cn) Qi^n) + a(en) Q(Cn " ??) = 0, (2.27) 

(Q(en),Q(^n-r?)) ^(0,0). (2.28) 

Similarly the right eigenstate with eigenvalue r(A) is a right separate state of the 
form (I2.22P for the same function Q(X). 

Remark 2.1. It is important to underline that, even if the system of equations (|2.27l) 
does not uniquely define the function (5(A), they nevertheless uniquely define the 
corresponding eigenstate up to a global normalization factor. 

Remark 2.2. The conditions (|2.26l) for the transfer matrix eigenvalue t(A) imply 
that a function Q{X) associated with r(A) through the system (I2.27p satisfies in fact 
the following larger systems of 2N equations: 

l^nh) QiCn hnlj) — ®(^n hnh) QiCn {hn 1)^) 

+ d(Cn - Kri) Qi^n - {hn “ 1)?/), Vu G {1,..., A^}, yhn G {0,1}. (2.29) 

Conversely it is easy to see that, if one can exhibit two functions r(A) and Q{X) 
satisfying (|2.29p and (|2.28p . then r(A) automatically satisfies the conditions (|2.26p . 
It means that the set of quadratic equations (I2.26p for t(A) in Theorem 12.11 can 
equivalently be replaced by a condition on the existence of a function Q{X) satisfying 

(ESD-dZSHI). 

2.3 Baxter functional T-Q eqnation 

The system of equations (12.291) which, due to Remark 12.21 leads to the character¬ 
ization of the transfer matrix spectrum and eigenvectors, is strongly reminiscent 
of Baxter’s famous T-Q equation [3]: it appears naturally as a discrete version, 
evaluated at the inhomogeneity and shifted inhomogeneity parameters only, of the 
following functional equation: 

r(A) Q(A) = -a(A) Q{X - g) + d{X) Q{X + p). (2.30) 

In fact if, for a given function t(A), there exists Q{X) solution of (I2.30p . then it 
is obvious that Q{X) satisfies (|2.29l) . The converse is of course not automatically 
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truqj. It would nonetheless be very convenient to be able to completely characterize 
the transfer matrix spectrum and eigenstates in terms of solutions of the continuous 
equation (|2.30p rather than of the discrete ones (|2.29p : the entireness condition 
for t(A) could then simply be rewritten in terms of Bethe-type equations for the 
roots Xj of the function (5(A), which would enable us to study the homogeneous and 
thermodynamic limit of our model in a rather standard way. 

Hence the whole problem, to pass from the discrete (|2.29p to the continuous 
(|2.3Up picture, is to understand whether there always exists, for each eigenvalue 
t(A) as characterized from the SoV approach, a solution Q{X) to the functional 
equation (|2.3Up which moreover satisfies (|2.28p . In addition, one has to be able 
to characterize the functional form of this solution, keeping in mind that the latter 
should be relatively independent from the corresponding eigenvalue t(A) if one wants 
to be able to reformulate our characterization of the spectrum and eigenstates in 
terms of a system of Bethe-type equations. 

Depending on the model one considers, this may be a difficult probleirH to pre¬ 
cisely identify the class of (5-solutions to (|2.30p associated with the transfer matrix 
eigenvalues t(A). However, in the present case, the situation is quite simple since 
we have to deal with polynomials of degree less than N only. One can therefore 
formulate the following theorem, which provides an alternative description of the 
transfer matrix spectrum with respect to Theorem 12.11 

Theorem 2.3. Let the inhomogeneity parameters satisfy the condition (12.lip . Then, 
a given function r(A) is an eigenvalue of the antiperiodic transfer matrix (|2.5li if and 
only if it satisfies the two following conditions: 

1. it is an entire function of X, 

2. there exists a polynomial Q of some degree R < N, 

R 

Q{X) = ll{X-Xa), (2.32) 

a=l 


for some set of roots Xi,...,Xr such that Xa / ft, Va G {!,...,i?}, V6 G 
{!,... ,iV}, satisfying with r(A) the functional equation (I2.30p . 

Whenever it exists, such a polynomial Q{X) is unique. 

Proof. It is quite straightforward to show that, if the conditions 1 and 2 are satisfied, 
then r(A) is an an eigenvalue of the antiperiodic transfer matrix: on the one hand 

^All that we can say in general is that Q(A) satisfies some inhomogeneous version of (12.3011 of 
the form 

r(A) Q(A) = -a(A) Q(A - y) + d(A) Q(A + y) + Fq{X), (2.31) 

with an additional term Fq{X) vanishing at all the points and — y, n € {1,..., A^}. 

^See for instance the recent papers IZ3EE1 which tackle this problem for the slightly more 
complicated XXZ and dynamical 6-vertex antiperiodic models: it appears that in these cases the 
functional form (and notably the quasi-periodicity properties) of the solutions Q{X) differ from those 
of t{X). The problem seems even more complicated when one considers spin chains with boundaries, 
so that it has been suggested in [96] to reformulate the SoV characterization of the transfer matrix 
spectrum in terms of solutions of a particular inhomogeneous functional equations of the type (IT3T1I 
rather than of the homogeneous one. 
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it is obvious from the previous discussion that the condition 2 of Theorem 12.11 is 
satisfied; on the other hand the fact that Q(A) is a polynomial of maximal degree N 
implies that —a(A) Q(A — r?) + d{X) Q{X) is a polynomial of maximal degree 2N — 1 
so that, from the entireness condition of r(A), the condition 1 of Theorem [2T] is also 
satisfied. 

Let us now show that, if the equation (I2.30p admits a polynomial solution Q{X) 
of the form ()2.32p for a given function r(A), then this solution is unique. Indeed, let 
us assume that there exists two different polynomial solutions P{X) and Q(A) to the 
equation ()2.30l) associated with the same function t(A). It means that 

-a(A)P(A - r?) + d(A)P(A + 77 ) -a(A)Q(A - r?) + d{X)Q{X + r?) 

P(A) ■ Q(A) 

so that 

-a(A) Wp,q{X) = d{X) Wp,QiX + ri), 

where Wp^q{X) stands for the quantum Wronskian of these two solutions: 

Wp,q{X) = Q{X) P{X - 77 ) - P(A) Q(A - 77 ). 

Taking into account that a(A) = d{X + 77 ), and using the fact that Wp^q{X) is a 
polynomial in A, we obtain that 

Wp,Q{X)=wp,Q{X)d{X) (2.36) 

where wp^q{X) is a polynomial in A which moreover has to satisfy the following 
quasi-periodicity condition: 


(2.33) 

(2.34) 

(2.35) 


wp^Q{X + r]) = -wp^q{X). 


(2.37) 


Since the only polynomial in A which is periodic of period 2rj is a constant, and since 
the only constant which satisfies wp^q = —wp^q is zero, one has that Wp^q{X) = 0. 
Hence Q{X) = P{X) once we have chosen the highest coefficients of both polynomials 
to be equal to 1. 

Let us finally prove that, for any eigenvalue r(A) of the antiperiodic transfer 
matrix (j2.5p . there exists a polynomial solution Q{X) of degree R < N and which 
does not vanish at the points ^j, 1 < j < N, to the functional equation (I2.3np 
associated with r(A). We recall that, if Q{X) is a polynomial of maximal degree 
N, then r(A) Q{X) and —a(A) Q(A — 77 ) -|- d{X) Q{X + 77 ) are both polynomials of 
maximal degree 2N — 1. Hence these two polynomials are equal if and only if they 
coincide for 2N different values of A, for instance at the 2N points and — i] for 
a G {1,..., A^}, i.e., if and only the following system of 2N equations is satisfied: 


r{Ca) QiCa) = -a{Ca) Q{^a “ r]), 
r{Ca - 77 ) QiCa - 77 ) = d{^a “ V) QiCa), 


VaG{l,...,W}. 


(2.38) 


Note at this point that, due to (I2.26p . the above system is in fact equivalent to the 
system of only N equations given by the first line of (I2.38p . Moreover, saying that 
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Q{X) is a polynomial of maximal degree N is equivalent to saying that Q{X) can be 
written in the following form: 


N+lN+l , _ 

= E IIfvT 

a=l b=l 
b^a 

In (|2.39l) . ^N+i is an arbitrary complex number, different from ... ,^n, which can 
be chosen at our convenience. Hence the system (I2.38p is equivalent to a homoge¬ 
neous linear system of N equations for the 2N -|- 1 unknowns Q(^i),..., (^(CAr+i), 
which can alternatively be thought of as an inhomogeneous linear system for the N 
nnknowns (5(^i), • • •, Q{Cn) in terms of the {N + l)-th one 

^ p -F - 

[Cr('^iV-Hl)]ab Qiib) = - n T- Q{iN+l)- (2.40) 



The elements of the matrix CT(^Ar+i) of this linear system are 

Vo, b G {1, 


r„ /t.. M _ .s: — — 


a(^a) 


C=1 

c^a 




,N}. 


(2.41) 


The determinant of this matrix is a rational function of ^jv+i which is not identically 
zero, hence it is possible to chose ^n+i for this determinant to be finite and non-zero. 
Then, for any given choice of Q(^Ar+i) 7 ^ 0, there exists one and only one nontrivial 
solution (< 5 (^ 1 ),..., Q(^n)) of the system (12.401) . which is given by Cramer’s rule: 


Q(^j) = Q(^n+i) 


with matrices Ct\^jv+i) defined as 


detAT [cl^^(gjv+i)] 
detAr[cT-(^Ar+i)] 




(2.42) 


N 


= (1 - Sbj)[Cr{CN+l)]ab “ h,j JJ 


i=l 


^N+1 - ’ 


(2.43) 


for all a,b & {1,..., N}. Note that the determinant of the matrices are also non¬ 
zero rational functions of ^at+i, so that it is also possible to fix ^Af-i-i such that 
none of them vanish. In that case one is ensnred that Q{^j) 7 ^ 0, Vj G {1,... ,1V}. 
Hence from (I2.39p we have obtained a polynomial of degree at most N satisfying the 
fnnctional equation (|2.30p and the reqnirement that its zeros do not coincide with any 
of the inhomogeneity parameters. It is important to remark that we have no criteria 
implying that the degree of Q(A) is exactly N: being described by an interpolation 
formnla in iV -|- 1 points, it can be a polynomial of any degree R < N. □ 

Hence, from Theorem 12.31 we know that there exists a bijection which relates 
the set of the transfer matrix eigenvalnes r(A) to the set of polynomials Q(A) of the 
form (I2.32P with maximal degree N: the image of t(A) is provided by the unique 
solution (5(A), in the aforementioned set of polynomials, to the functional equation 
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(|2.30l) associated with t(A). From now on we shall denote this unique solution by 
Qr{X), and we shall denote the corresponding right and left eigenstate by | Qr ) and 
(Qr I respectively. 

The Baxter equation as usual leads to the Bethe equations for the roots of the 
polynomial Q{X): 


O'iXa) TT- Xa — Xb — V _ 2 
d{Xa) Xa — Xb + r] 


Va E {1,... ,i?}, 1 < R< N. 


(2.44) 


Theorem l2. 31 therefore ensures us that the solutions to these Bethe equations provide 
a complete description of the spectrum (and eigenstates) of the antiperiodic transfer 
matrix (ESD. 

Before concluding this subsection, we would like to mention that it is possible to 
obtain an alternative (but equivalent) description of this spectrum, as stated in the 
following theorem. 


Theorem 2.4. Let M = N/2 if N is even or M = {N — l)/2 if N is odd. Then 
t(A) is an eigenvalue of the antiperiodic transfer matrix (12.51) if and only if it can 
be written in the following form, 

'r(A) = [p(A - rj) q{X + r])- q{X - rj) p{X + rj )], (2.45) 

with q{X) a polynomial of degree R < M and p{X) a polynomial of degree N — R such 
that 

\ b(A) q{X - 7?) + q{X) p{X - rj)] = d{X). (2.46) 

Proof. Let r(A) be an eigenvalue of the transfer matrix. Then Theorem 12.31 implies 
that there exists a (unique) polynomial Qt{X) satisfying the Baxter equation (12.30^ 
together with t(A). From Theorem 12.11 it is easy to see that —t(A) is another 
eigenvalue of the transfer matrix, so that there also exists a polynomial Q-r{X) 
which satisfies (j2.3np together with —r(A). We then define q{X) € {Qr(A), Q-r(A)} 
to be the polynomial with the smaller degree and p{X) € {Qt(A), Q_r(A)} to be the 
other one. If the two polynomials Qt{X) and Q-r{X) have the same degree, we fix 
for instance q{X) = Qt{X) and p{X) = Q-r{X). 

Using the fact that q{X) and p{X) satisfy the Baxter equation with opposite 
eigenfunctions, we obtain the identity 

-a(A) q{X - r/) + d{X) q{X + tj) _ a(A) p{X - rj) - d{X) p{X + rj) 

q{X) p{X) ’ ^ ^ 

or equivalently, 

a(A) W,,p{X) = d{X) W,,p{X + r?), (2.48) 

where we have defined the function Wq^p{X) as 

^ b(A) qiX -v) + qW P(A - r?)] • (2.49) 


Since Wq^p{X) is a polynomial in A and since a(A) = d{X + p), the equation (|2.48p is 
satisfied if and only if 


Wp,q{X) = d{X). 


(2.50) 
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It also means that the degree of q{X) is less or equal to M. 

Let us now show that the transfer matrix eigenvalue t(A) is of the form (I2.45p . 
By definition we know that 

r(A) q{X) = e[-a(A) q{X - rj) + d{X) q{X + rj)] (2.51) 

where e = 1 if q{X) = Qr{X) and e = —1 if q{X) = Q-riX). Using now that 
a(A) = Wp^q{X + rj) and that d{X) = Wp^q{X), we obtain 

r(A) q{X) = e [g(A) p{X - rj) q{X + rj) - q{X) p{X + rj) q{X - rj )], (2.52) 

which implies (12.451) . 

Vice versa, let q{X) and p{X) be two polynomials of degree R and N — R respec¬ 
tively and which satisfy (|2.46p . Let r^'''^(A) and t^~\X) be equal to the right hand 
side of (j2.45l] with -|- or — sign respectively. We can first remark that, from the 
definition (j2.45p . r(^^(A) are obviously polynomials in A of degree N — 1. Moreover, 

r(=^)(A) q{X) = ± ^ [p(A - rj) q{X + rj) - q{X - rj)p{X + rj)] q{X) 

= I ^ “ 'd) qW + pW 9(A - rj)] q{X + rj) 

- ^ [p{X + rj) q{X) + q{X + rj) p(A)] q{X - r?) | 

= ± [diX) q{X + rj)- o(A) q{X - rj)]. (2.53) 

Similarly one can show that 

(A) p(X) = =F [ - a(A) p(X - rj) - d(X) p(X + p)] ■ (2.54) 

Hence t^~^\X) = —)(A) satisfies the functional equation (12.30^ with q(X) whereas 
-1(A) = —T^~^\X) satisfies the functional equation ()2.30p with p(X). Finally it 
is easy to see that, given any inhomogeneity parameter ^j, p(X) and q{X) cannot 
both vanish in — p: this would imply from (|2.46l) that d(^j — p) =0, which is 
obviously not true. Hence, taking also into account (12.531) and (12.541) . this means 
that T^~^\X) and t^~\X) both satisfy (j2.26p . so that they are both transfer matrix 
eigenvalues. □ 

2.4 ABA-type representations for the transfer matrix eigenvectors 

In this section we present an alternative way to write the separate states, and hence 
the eigenstates of the antiperiodic transfer matrix, in a form which is strongly remi¬ 
niscent of the form of the Bethe states as obtained in the framework of the algebraic 
Bethe ansatz. It is worth remarking that this type of rewriting can be in fact derived 
for a large class of integrable quantum models solvable by SoV method. 

Let us first define a simple separate state, the state (11 associated with the 
constant function q;(A) = 1. The corresponding right separate state | 1) is defined 
similarly using (|2.22p . Then, if a(A) is a polynomial, the separate state (a | or | a) 
can be obtained by multiple action of the “creation” operator D, evaluated at the 
roots Uj of a, on the separate states (11 or | 1). 
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Proposition 2.1. Let a{X) be a polynomial of the form, 

R 

a(A) = [](A-afc). (2.55) 

k=l 

Then the corresponding separate states {a \ and \ a) can he written as 

R R 

(a| = (-l)^'^(l|[]ll(afc), |a) = (-1)^"^ J]Zl(afc)|l). (2.56) 

k=l k=l 

Proof. The proof is straightforward using (I2.16P and ()2.17p . □ 

The representation (|2.56l) can be used notably for the eigenstates | Qr ) and ( \ 

corresponding to the eigenvalue r(A) of the transfer matrix, by means of the unique 
polynomial Qt(A) which, in virtue of Theorem 12.31 solves the functional equation 
associated with t(A). Since the roots of Qt(A) can be obtained as the solutions to 
the corresponding Bethe-type equations, the analogy with algebraic Bethe ansatz 
is then particularly obvious. It is also possible to represent these eigenstates in a 
slightly different form. Let us to this aim define some other simple left and right 
separate state that we shall denote by (lait | and | lait) respectively, and which are 
given by a function lait (A) which alternates sign between the inhomogeneity and 
shifted inhomogeneity parameters, namely 

lalt(ea) = l, lalt(ea - r?) = -1, Vo E {1, . . . IV}. (2.57) 

This second representation for the eigenstates | Qr ) and ( Qr \ corresponding to the 
eigenvalue r(A) of the transfer matrix then uses the polynomial Q_t-(A) which solves 
the functional equation (j2.3np for the eigenvalue —t(A) of the transfer matrix: 

Proposition 2.2. Let r(A) be an eigenvalue of the antiperiodic transfer matrix (12.5p 
and let 

R N-R 

Qr{X) = lliX-Xk), Q-r{X)= HiX-Xk), (2.58) 

k=l k=l 

be the unique solutions to the functional T-Q equation (I2.30p associated with the 
eigenvalues r(A) and —t(X) respectively. Then the transfer matrix eigenstate \ Qr ) 
with eigenvalue t(A) can he represented in following forms: 

R 

|q^) = (_1 )«a^]JZ)(A,)|1 ) (2.59) 

k=l 

R 

n d{^a) N-R 

= ^_i)iN-R)N k=i^ I (2.60) 

n d(Xa) ^=1 
k=l 
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3 Scalar products of separate states 


In this section, we explain how to compute the scalar products of separate states of 
the form (j2.2ip and (|2.22p . In general, the latter can by construction be represented 
in terms of the determinant of a weighted sum of two Vandermonde matrices in¬ 
volving the inhomogeneity parameters of the model. We shall notably focus on the 
case where one of the two separate states is an eigenstate of the antiperiodic transfer 
matrix. We shall see that in this case it is possible to represent the corresponding 
scalar product in a more convenient form for the study of the homogeneous and 
thermodynamic limit, namely, similarly as what happens in the ABA framework 
for the scalar product of an on-shell and an off-shell Bethe states, in terms of the 
determinant introduced in |83] (Slavnov determinant). 


3.1 General determinant representation for the scalar product of 
two separate states 

Let us start by recalling the general determinant representation for the scalar prod¬ 
ucts of two separate states. Note that this type of representation is a direct con¬ 
sequence of the factorized form of these states in the SoV basis, and hence can be 
shown for a large variety of models solvable by SoV. 


Theorem 3.1. Let (a | be a left separate state of the form (12.211) . and let \ /3) be a 
right separate state of the form (|2.22p . Their scalar product can be written as 


{a\P) 


detAT 

^(Ul) 


(3.1) 


where corresponds to the following weighted sum of two Vandermonde ma¬ 

trices: 


= ^ Ol{ia) ^{ia) + {fa - hf ^ 0:{fa - V) Pifa - V) 

a{fa)a{fa -V)l3{fa -v) 


= a{^a) ma 


cb-l_ 
Sa 


d{fa - h)«(^a)/3(?a) 


iCa - V) 


b-1 


(3.2) 

(3.3) 


The proof of this theorem is straightforward. It was given in |46] in the case of 
a more general XXZ model. 

Note that it is possible to use this determinant representation to show the (ex¬ 
pected) orthogonality of the eigenstates of the antiperiodic transfer matrix, which 
therefore form an orthogonal basis of the quantum space of states of the model. 


Corollary 3.1. Let t(A) and t'{X) be two different eigenvalues of the antiperiodic 
transfer matrix (12.5p . Then the corresponding eigenstates are orthogonal: 


{Qt' \ Qt) = {Qt \ Qt' ) = 0 . 


(3.4) 


Proof. r(A) and t'(A) are two different polynomials of degree iV — 1: 

N N 

= t'{X) = (3.5) 

6=1 6=1 
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It means that the vector with components Vf, = Cf, — is non trivial, and it is easy 
to check that 

N 

= (3.6) 

b=l 

□ 

The aim of this section is to show that it is possible to rewrite the expression (EH) 
for the scalar product of two separate states into some more convenient forms for the 
consideration of the homogeneous and thermodynamic limit. We shall in particular 
link this formula to some (generalizations of some) other determinant representations 
that have already appeared in the literature in the framework of algebraic Bethe 
ansatz; the Izergin determinant, first introduced as a representation for the partition 
function of the six-vertex model with domain wall boundary conditions [82] . and 
the Slavnov determinant, which appears in the ABA framework as a convenient 
representation for the scalar product of an on-shell and an off-shell Bethe vectors 
[83] . To this aim, we shall first derive some identities that will be used for these 
reformulations. 


3.2 Some useful identities 

As announced above, we want to transform the representation EH into a more 
convenient form for our purpose. To this aim we introduce, following [21], some 
convenient notations that we shall use throughout all this section. For any set of 
complex numbers {x} = {xi,..., xm}, and a function /, we define 


N 


= n 

n=l 


y -Xn±r] 
y-Xn 



detM [xl ^ - f{Xa) {Xa ± yf 

VM) 


(3.7) 

(3.8) 


Note that the representation (I3.ip - ([3.3I) of the scalar product can easily be rewritten 
by means of the notations (13.71) and (|3.8I) : it is enough to choose the functions a (A) 
and /3(A) to be some polynomials (which by interpolation is always possible) that 
we express in terms of their roots {ai,..., an} and {/3i,..., Ps} as 


R s 

«(A) = n(A-a„), /3(A)= n(A-/?m), (3.9) 

n=l m=l 


to obtain that 


N 

(«i/ 3 )=n A 


{?} 


n=l 




(3.10) 


We shall reformulate this expression by means of a few identities, involving the 
quantities (EH and (13.81) and their relations to the Izergin and Slavnov determinants, 
and that we now derive. 
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We start with some preliminary identity which relates the two functionals j [/] 
and (|3.8I) when /' is related to / via a ratio of the two functions and 

E,-, EH: 

Identity 0. For any set {x} = {xi ,..., xm} of arbitrary complex numbers, we have 


■^{x} [■/’] “ -^{x} 


E 


{x} 




/ 


(3.11) 


Proof. To prove this identity we use an expression for the determinant of the sum 
of two matrices, 


V{{x}) ^ 


{x} 


E+ 

1^1 X 

Efx 

L J 


1 M 


M 


E'" E n -/(^-)n 


m=l 


Xn-Xm + rj 


V{{x — hr]}) 


hi=0 hM=0ri=l \ 

1 1 / M \ 

= E ■ ■ ■ E ( n ) ^({^ + ^h}) 

hi=0 hM=0 Vn=l / 

= F{{,U).4+,[/]. (3.12) 

We have also used here the identity (j2.23p for the Vandermonde determinants. □ 


We shall now formulate a generalization of a result of Kostov [8^IWl[9nllMl[W| 
concerning the relation between the functionals [/] and the Izergin determinant 
(see also [93]). 

For ]j. G C and two sets {x} = {xi,...,XAr} and {y} = {yi,... ,yN} of arbi¬ 
trary complex numbers, we introduce the following function, that we call generalized 
Izergin determinant: 


N 


i!i‘'(w.{s}) = 


n (Xa -Vb + V) 

del [t^(xa - ?/b)], 


-N V{xi,...,XN)V{yN,---,yi) “n 

where we have defined 


tf,{x) = -- 


1 


(3.13) 


(3.14) 


X X + T] 

For /i = 1 this formula gives the partition function of the rational six-vertex model 
with domain wall boundary conditions 


Identity 1. For ]j, G C and any sets {xi,...,XAr} and {yi,... ^y^} of arbitrary 
complex numbers, we have 




{y} 


= i-^rK 


{y} 


]iE_ 


{x} 


(3.15) 

(3.16) 
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Proof. We shall prove here the first equality, the second one can be proven in the 
same way (or by replacing rj by —rj in the Izergin determinant). 

Let us first introduce the following auxiliary matrix with elements, 

N 

Mab=Y\_{Xa-yk + V)ti,{Xa-yb), (3.17) 

k=l 


SO as to rewrite (j3.13p as 


detw 7V4 


V{xi,...,XN)V{yN,---,yi)' 

We also introduce a set of N polynomials of degree N — 1, 

N 


Zb{x) = Y\{x - yk + T]) = '^Cj^bX^ S 


l<b< N. 


(3.18) 


(3.19) 


k^b j=l 

The coefficients Cj^b of these polynomials can be seen as the entries of an x 
matrix C. It is easy to observe that 




y,Zb{Xa 


N 


^)n 


Xg-yk + y 
Xg - yk 


Zb{Xg) 



(3.20) 


which means that Ad can be factorized into a product of two matrices, so that 


det Ad 

N 


(-1)" r({i)) 



detC. 

N 


(3.21) 


It now remains to compute the determinant of the matrix C. It is easy to see from the 
definition of the coefficients Cj^b that the product of C with a Vandermonde matrix 
is diagonal: 

N 

= KbYliya - yj), (3.22) 

i=i j¥=o. 

which means that detwC = V{yN, ■ ■ ■ ,yi)- □ 


Note that the Izergin determinant (I3.13P is defined only in the case where the 
cardinality of the two sets of parameters {x} and {y} are equal, so that Identity [U is 
a priori valid only in this case. To overcome this restriction, it is however possible 
to generalize the part of Identity [T] which concerns the relation between the two 
functions and All'll to cases in which the number of parameters 

X and y does not coincide. 

Identity [Ij. For two sets {xi,..., xm} o,nd {yi,..., yw} of arbitrary complex num¬ 
bers and y / 1, we have 




{i-y) 


N-M^- 


{x} 




{?/} 


(3.23) 
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Proof. The case N = M follows immediately from Identity [TJ Let us therefore 
suppose that M > N (the case M < N can be proven in the same way startiM from 
the right hand side). The proof is based on the following simple observatioiu, 


A + 

•^{y} 




{x} 


= (1 -/^) 


N-M 


lim 


At 




fiE- 


from which we can use Identity [T] and obtain 


A 


{y} 




N-M 


lim Af 1 




{y}U{vi,...,VM-N} 


Computing the limits we obtain (13. 2, ip . 


(3.24) 

(3.25) 

□ 


Note that for N > M both sides of (j3.23D are polynomials in ^u, so that the 
identity holds also in the case ^ = 1. This leads to the following simple observation: 


Corollary 3.2. For arbitrary sets {xi,..., xm} and {yi,..., of complex num¬ 
bers such that N > M, we have 




= 0 . 


(3.26) 


The next identity that we shall derive will enable us to relate the scalar product 
of two separate states, one of them being an eigenstate of the transfer matrix, with 
the Slavnov formula [83] for the scalar product of an off-shell and an on-shell Bethe 
vectors in the framework of the algebraic Bethe ansatz. 

For /r G C, for two sets of parameters {x} = {xi,..., xm} and {y} = {yi ,..., pm}, 
and a set of inhomogeneity parameters {^} = {^i,... ,^ 7 v} (which may possibly co¬ 
incide) with M < N, we introduce the following function, that we call Slavnov 
determinant: 

M 

n (xj -Vk+r]) 

4?’(U).WI{«}) = (3.27) 

V{xi,...,XM)V{yM,---,yi) M 

In (13.271) . the matrix is defined by its elements as 

[7^i'"i({x},{y}|{C})]^.^ = ^E^^.^{yk)t{xj - Vk) - ^=^-^t{yk - Xj), (3.28) 
where t{x) = ti(x) is given by (|3.14l) . 

As already mentioned, this formula appeared initially in the ABA framework as a 
representation for the scalar product of a Bethe state associated with generic param¬ 
eters {y} (off-shell Bethe state) and a Bethe state associated with Bethe roots {x} 
satisfying a system of Bethe equations (on-shell Bethe state) for an inhomogeneous 
model with inhomogeneities ^i,..., and with a twist p: 

Et 1 (Xm) 

I^Eyx„.) = (3.29) 

Note at this point that the Bethe equations (|2.44l) introduced in the previous section 
correspond to a special case of (I3.29D with p = —1. 

^Here and in the following we mean that the limits are taken independently one at the time. 
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Identity 2. For fj, G C, for {xi,... ,xm} a solution to the Bethe equations (|3.29l) 
associated with a set of inhomogeneity parameters {^i,..., f,N}, and for {yi,..., j/m} 
a set of arbitrary complex numbers, we have 


{!/}!{«})=-^wub) 



(3.30) 


Proof. Let us first suppose that the parameters Xj are pairwise distinct, and let Q{\) 
be the normalized polynomial of degree M with roots xi,, xm- 


M 


Q(A) = n(A-x,). 

i=i 


(3.31) 


We define from Q{X) M different polynomials Qk{X) of degree M — 1 as follows: 


Qfc(A) = ]J(A-Xj), l<k<M. (3.32) 

Then, similarly as what has been done in the proof of Identity [H we can define a 
2M X 2M matrix C from the coefficients of the following 2M polynomials: 

2M 

Qfc(A)gfc(A + 7?) = ^Cfc,aA“-\ l<k<M, (3.33) 

a=l 

2M 

Qfe(A)Q(A + y) = ^Cfe+M,aA“-\ l<k<M. (3.34) 

a=l 

Evidently the coefficients Ck, 2 M all vanish for k = 1,..., M, but the 2M x 2M matrix 
C is invertible, and it is not difficult to compute its determinant using the following 
identities: 

2M 

^ ^ Cb,a (^fc h) ^b,k Qki^^k) y), (3.35) 

a=l 

2M 

^ ^ Cb,a ^b,k+M Qi^^k T ^) Qfc(®fc) T ^b,k QkiFk) Qki^k T h)- (3.36) 

a=l 

It means that the product of the matrix C with the Vandermonde matrix constructed 
from the variables xi — y,..., xm — j ■ ■ ■, xm is triangular, and therefore the 
determinant of the matrix C is given as 


detC 

2M 


M 

n [Qi (Xk) Qk{xk - r]) Q{xk + r])] 

k=i _ 

V(xi-r],...,XM-r],xi..., xm) 


M 

AT(A/+1) ^ ^ |- -| 

= (-1) 2 II [Qk{xk) Qk{xk + ^?)J • 

k=l 

On the other hand let us compute the following product. 


A 


{a:}U{y} 



detC 

2M 


^ ( gAP g{i,2) X 

detsM ^(2,1) ^(2,2) j 
V{{x}U{y}) 


(3.37) 


(3.38) 
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where the matrix in the numerator is presented in a block form, {1 < p,q < 2) 

being M x M matrices which can easily be computed. For instance, the upper left 
block ];ias for elements 

2M 

- f^E+^ixk) (Xk - vT~^) 

a=l 

= Qj{xk) Qj{xk + v) - iJ^E^^^ixk) Qj{xk) Qj{xk - i) 

= Sjk Qkixk) (Qkixk + v)- ^J'E^^-^{xk) Qkixk - ?/)) , ( 3 . 39 ) 

so that ^dd) = Q (iue Bethe equations (|3.29D . It means in particular that we 

do not need to compute the lower right block have access to the quantity 

(|3.38l) . The lower left block is a diagonal matrix, 

2M 

(xl-^ - /rF;+j(Xfc) (Xfc - vT~^) 

a=l 

= Qjixk) Q{xk + ri)- yLE^^^{xk) Q{xk) Qj{xk - r]) 

= SjkQk{xk)Q{xk + ri)- (3.40) 

Finally, the upper right block jg given by 

2M 

{yt~^ - fkE+^ivk) ivk - 

a=l 

= Qjivk) Qjivk + 'n)- Qjivk) Qjivk - v) 

= -^Q{yk)Q{yk-v) ['H^''^({x},{y}|{^})]^.^. (3.4i) 

It means that 


M 


A 


{x}U{y} 


yEt 


{«} 


n [Q{yk) Qivk - y) Qk{xk) Qk{xk + y)] 


detC = 
J 2M 


k=l 


y({x}U{y}) 


X det^({x},{?/}|{C}), (3.42) 

M 


which leads directly to the identity (j3.30p . 

Note that both 5^^({x}, {?/}|{^}) and formally contain zero 

over zero terms if computed for two or more coinciding Xj. So that the identity 
(|3.27l) has to be meant as a limit to the Bethe roots in the case of two or several 
coinciding roots. In fact we can introduce a set of parameters {xd)} pairwise distinct 
for e > 0 and converging to the solution of the Bethe equations {x} when e —?• 0. 
We have to prove that the l.h.s. and r.h.s. of (|3.27p are finite and coincide under 

3.231 being proven for arbitrary values of the parameters. 


this limit. The identity 


implies that A 




yE. 


{?} 


has a smooth limit for 


0 . 


Now we can repeat 

the proof above developed for pairwise distinct roots checking that at any step all 
remains finite and that one just reproduces 5^^({xd)}, {?/}|{^}) for e ^ 0. □ 
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If the number M of parameters x is exactly half of the length of the chain we 
can directly apply Identity [1] and Identity [2] to obtain the following result: 


Corollary 3.3. Let N = 2M. Then, for /r G C, for {xi,... ,xm} « solution 
to the Bethe equations (I3.29P associated with a set of inhomogeneity parameters 
{?!) ■ ■ ■) ^n}, and for {yi,..., yM} a set of arbitrary complex numbers, we have 









(3.43) 

(3.44) 


This relation is of prime importance for the computation of the scalar products 
of separate states. The inconvenient restriction on the number of parameters x and 
y can be relaxed thanks to Identity m which enables us to formulate the following 
generalization of (j3.43p : 


Corollary 3.4. Let {xi,... ,xm} be a solution of the Bethe equations ()3.29p asso¬ 
ciated with the set of inhomogeneity parameters and the twist y, and 

let {yi,... ,yM} be a set of arbitrary complex numbers. Then the following relation 
holds: 


5i^({x},{?/}|{e}) 


(1 - 


_l^^{x}U{y} 


(3.45) 


Finally, we would like to generalize the previous results to cases for which the 
number of parameters x and y are not obligatory the same. To this aim we introduce 
the following generalization of the Slavnov determinant (|3.27l) . which is defined for 
two sets of parameters {x} = {xi,..., xm} and {y} = {yi ,..., yu-es}, with 5* > 0, 
a set of inhomogeneity parameters {^} = ..., and a twist y: 




M M-es 

n n (xj -vk+h) 

j=i k=i _ 

V{xi,..., xm) V (yM-hS, • • •, yi) 

X det ^('")({x},{y}|{C}), (3.46) 


with 

[H(^)({x},{?/}|{C})] = yE+(yk)t(xj - yk) - } t{yk - Xj), if j < M, 

^{x}yyk) 

= yE-^^pk) ivk + if j > M. 

^{xyyyk) 

This type of object was first introduced in [92]. This enables us to formulate the 
last identity of this subsection: 


Identity 3. Let {xi,..., xm} be a solution of the Bethe equations ()3.29p associ¬ 
ated with the set of inhomogeneity parameters {^i,... and the twist y, and let 
{yi, ..., yu-es} be a set of arbitrary complex numbers. Then, 




A 


{x}U{y} 



(3.47) 


The proof of this identity follows the same lines as the proof of (|3.30l) . It is 
however more cumbersome, so that we give its details in Appendix A. 
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3.3 An alternative representation for the scalar product of two sep¬ 
arate states 

We shall now use the identities that have been derived in the previous subsection to 
rewrite the representation (13. ip (or (13.1011 ) for the scalar product of two generic sep¬ 
arate states (a \ and | /3) into a form for which the consideration of the homogeneous 
limit is completely straightforward. 

Let us first notice that Identity [0] enables us to rewrite (j3.inil into a slightly 
simpler form: 

Proposition 3.1. Let a{\) and /3(A) be two polynomials of respective degree R and 
S, and which can he factorized as 

R S 

a{X) = [] (A - a„), /3(A) = [] (A - /3™), (3.48) 

n=l m=l 

in terms of some sets of roots {ai,... ,aR} and {/3i,... ,/35}. Then the scalar product 
of the corresponding separate states {a \ and \ [3) can be written in the following form: 


{a\P) = (-l)^(«+5) d{aj) JJ d{/3k) 

j=i k=i 


-E 




, (3.49) 


where we have used the notations (isip-dsSD. 

We can now use the identities [T] and m and we obtain the following result: 


Theorem 3.2. Under the same hypothesis and notations as in Provosition AS.R the 
scalar product of the two separate states {a \ and \ j3) can he written as: 


R S 

{a\/3) = (-1)^(«+^) J] d{a,) d(/3fc) 

j=l k=l 


X A 




-E 


{«} 


. (3.50) 


In the particular case where R + S = N, the scalar product of the two separate states 
{a I and \ /3) can be written in terms of the generalized Izergin determinant (j3.13p 
as: 


R 


{a\(3) = (_i)^(-R+s+i) di/3k) 


7 = 1 


k=l 


^^({ai,---,ai?}U{/3i,...,/35},{^}). (3.51) 


Note that, as announced before, the representation (I3.50p is completely smooth 
with respect to the homogeneous limit in which all the inhomogeneity parameters 
fj, ^ < j < 3V, tend to the same value. 
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3.4 The scalar product of a generic separate state with an eigen¬ 
state of the transfer matrix 


We shall now focus on the case in which one of the two separate states is an eigenstate 
I Qr ) of the transfer matrix. As we shall see, the knowledge of such scalars products 
is indeed sufficient to obtain some adequate representations for the form factors of 
local operators. This fact has to be put in relation with what happens in the ABA 
framework, where the consideration of the scalar products of off-shell and on-shell 
Bethe vectors is used to compute the local spin form factors m- 

In the whole subsection, (a \ will denote a separate state associated with a given 
polynomial Q!(A) of degree M, 

M 

a(A) = ]J(A-a,), (3.52) 

n=l 

and I Qr) will denote a given eigenstate of the antiperiodic transfer matrix, associ¬ 
ated with an eigenvalue r(A). We recall that such an eigenstate can be written in 
two different forms, either by using the polynomial 

R 

Qr{X) = ll{X-Xk) (3.53) 

k=l 


satisfying the functional T-Q equation (l2.3Up with r(A) itself, or by using the poly¬ 
nomial 

N-R 

Q.r{X)= HiX-Xk) (3.54) 

k=l 

satisfying the T-Q-equation with —r(A) (see Proposition [221) • We also recall that in 
this framework the set of roots {Ai,..., Xr} of Qt(A) satisfies the Bethe equations 
()3.29l) with fi = —1 (as well as the set of roots {Ai,..., Xn-r} of Q-r{X)). 


Proposition 3.2. 

following forms: 

( « I Qr ) = ( 
= ( 


The scalar product (a | Q^) can be written in the two possible 


M R 

-l)N(R+M) d{an) n d{Xk) 

n=l k=l 

M N-R 

-lf(N-R+M) Y[d{an) n d(Xk)A 
n=l k=l 


-E 


{A}U{a} 


{«} 


E 


{A}U{a}J 


(3.55) 

(3.56) 


Proof. The first line is a direct consequence of (|3.49p . The second one can be obtain 
in the same way from the fact that 

QriCa ~ v) _ Q—rifa ~ V) /„ 

QriCa) ~ Q-rifa) ' ^ ^ 

□ 


We can now use the identities obtained in the previous subsection to prove the 
following results for the scalar product of the separate state (a \ with the eigenstate 

IQr). 
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Theorem 3.3. If M < R the scalar product between {a \ and \ Qr ) vanishes: 


{a\Qr) = 0. 


(3.58) 


If M = R the scalar product between {a \ and \ Qr) can he written either as an 
Izergin determinant (13.131) . 


M 


N-M 


(«iQr)=(n^K) n 4 '^(wu{A},{e}), (3.59) 


\n=l k=l 

or as a Slavnov determinant (13.271) . 

/ M 


(a|Q.) = (-!)“ 2 


M r,N-2M 


]Jd(an)d(An) (3-60) 


\n=l 


If M > R the scalar product between {a \ and \ Qr ) can be written as a generalized 
Slavnov determinant (|3.46p .' 


M 


R 


{a\Qr) = {- 1)^2 


R r,N-M-R 


n<*(“»)I (3-61) 


^n=l k=l 


Proof. The proof of this theorem is straightforward. To prove (|3.58p it is enough to 
use (j3.56p and Corollary 13.21 (I3.59p follows from (j3.56l) and Identity [H The two last 
representations follow directly from (j3.55D , and from the identities ^1 El and EJ □ 


It is important to mention that all the representations of Theorem 13.31 for the 
scalar product of the eigenstate | Qr ) with the generic separate state (a \ remain 
finite (and manageable) in the homogeneous limit. 

Note finally to conclude this section that the fact that we can represent the scalar 
product (a I Qr) as a Slavnov determinant allows us to take, as usual, the limit in 
which the two states are equal. We therefore obtain, as in ABA, a representation 
for the “square of the norm” of the eigenstate | Qr) (he. for the scalar product 
(Qr I Qt )) in terms of a Gaudin determinant. 


Corollary 3.5. The “square of the norm” of the eigenstate \ Qr) of the antiperiodic 
transfer matrix is given by 


R 


/ R n (Am ~ An + f?) 

{<3.1 <3.) = (n 

m^n 


R 


(3.62) 


with 


[<hr 


I m,n 




dXn I d(^\m) , ^ Xyyi Xfy + 7] 

\ 0=1 / 


(3.63) 
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4 Form factors of local operators 


In this section we compute the form factors, i.e. the matrix elements of the local 
spin operators between the eigenstates of the transfer matrix. 

Let us first consider the local operator a~. Its matrix elements between two 
eigenstates (Qr \ and | ) of the transfer matrix can be computed by acting with 

this operator on one of these two states, which can be done as usual by means of the 
solution of the quantum inverse problem [571165] . The reconstruction formula [461157] 
for this operator takes the form. 


a 


n 


{-If 




mn) (A m) \ 

\j=n+l / 


It enables us to formulate the following result: 


(4.1) 


Theorem 4.1. Let \ Qr) and \ Qr') he two eigenstates of the transfer matrix, with 


R 


Qr{X) = ll{X-Xk), 

k=l 


R' 


k=l 


(4.2) 


Then the corresponding form factors for the operator are given by the following 
expressions: 

• if \R — R'\ > 1, the form factor {Qr\a'~ \ Qr' ) vanishes, 

{Qr\a-\Qr') = ^-, (4.3) 


if R = R' + 1, the form factor {Qr\cr^ \ Qr') can be expressed in terms of the 
determinant of a matrix of size R as 


{Qr\a-\Qr') = 2 ^^-^^{-l) 


N-2R I 1^Ar-l Qriin) 


where 


and 


Qt' (^n) 


R R^ 

n «n(Afc) n dn{X'j) 
k=l j=l 

V{Xi,...,XR)V{f,,...,X[) ^R 


- det.F-({A},{A'},^„), (4.4) 


n N 

«n(A) = n(A-^,+r/) n 

j=l j=n+l 

n N 

dn{X) = ll{X-^j) n (A-e,+r?), 

j=l j=n+l 


(4.5) 

(4.6) 


a(A‘ 


^j,k ~ fxf Qt{X'i: — v) t{Xj — A).) + Qr{X'k + rj) t{X'i^ — Xj), for k < R, 


^j,R ~ ^(Aj 


(4.7) 
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if R' = R + 1, one has similarly 


( Qr I O-n I Qr' ) = 2 


N-2R' 


(- 1 ) 


N—1 Qr'iCn h) 


Qrifn V) 


R R' 

n an{\k) n 

fc=l_j=l 


detR-{{X'},{X},Cn); (4.8) 


finally, if R = R', the matrix element { \ \ Q^-/) takes the form 


( Qr I O-n I Qr' ) = 2 


R R' 

n an(Afc) n dn{X'j. 

N-2R-1 i -,\N-R Qr^fn) k=l j=l 


i-iy 


Qr'iCn) V{Xi,...Xr)V{X'j^,,...X[) 


X det 
R L 


.^({A}.{V}) + P<”I({A'),{A}.{„)1, (4.9) 


where the elements of the matrix J^({A}, {A'}) are 

^ 'A = Qri^'k ~ v) t{^j — ^'k) + Qri^'k + v) t{^'k ~ (4.10) 

and where ^n) is a rank one matrix, 

'P^k — I Qri^k — h) + Qr(Afc + r/)^ t{Xj — ^n). (4-11) 

Proof We consider the case R> R' (the case with R < R' can be done in a similar 
way). Using the reconstruction formula (|4.ip and the fact that {Qt\ and \Qt’) 
are eigenstates of the transfer matrix with respective eigenvalues t(A) and we 

immediately get the following representation for the matrix element (Qr \ erfi \ Qr )■ 

n—1 N 

n r{Cj) n r'iCj) 

{Qr\afi\Qr') = i-lf — -- ( Qr \ D{^n) \ Qr' )■ (4.12) 

n a(0) 

i=i 

We can now use the ABA type representations (I2.56p from which we straightfor¬ 
wardly obtain that 

D{in)\Qr') = {-lf\a), (4.13) 

with a;(A) being a polynomial of degree R' + 1 constructed from Qr'{X) as 

a{X) = {X-in)Qr'{X)- (4.14) 

Hence, it reduces the computation of the form factor ( Qt | o'n I Qr’) to the compu¬ 
tation of the scalar product {Qr\oi) between the eigenstate (Qr \ and the separate 
state I a ). We can therefore directly use the results of Theorem 13.31 More precisely, 
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in the case R > R' + 1, we can use the result for the scalar product (13.581) to show 
that the corresponding form factor is zero. If i? = i?' + 1 we apply (|3.6ni) which 
leads to the result (03]). Finally, in the case R = R', we can rewrite the gener¬ 
alised Slavnov determinant (I3.6ip as a determinant of a sum of two R x R matrices 

031 ). □ 

We can now use the symmetries dMl) of the transfer matrix to compute the 
matrix elements of the local operators and from the ones of a~. 

Theorem 4.2. Let \ Qr) and \ Qr') he two eigenstates of the transfer matrix con- 
strueted from polynomials Qt{^) and Q^i{X) with respective degree R and R'. The 
matrix elements ( Qr | <7^ | Qt' ) and ( Qr | o'n I Qr') of and are given in terms 
o/ ( Qr I O-n I Qt' ) aS 


( Qr I < I Qr' ) = 2 {R! - i?) ( Qr k" I Qr' ), ( 4 . 15 ) 

( Qr k+ I Qr' ) = (-1)^-^' ( Qr k" I Qr' ). ( 4 . 16 ) 

Proof. The statement of the theorem follows from the fact that we can simply obtain 
(T^ and (T+ from a~ by means of the operators and (1231): 

< = at = r^a-rP ( 4 . 17 ) 

Hence, to obtain the matrix elements of af and in terms of those of a~, we 
just have to compute the action of the operators and F^ on the corresponding 
eigenstates. Using the following limit 

we obtain 

S^\Qr) = {N - 2 R)\Qr). ( 4 . 19 ) 

The action of the operator F^ on the eigenstates of the transfer matrix can be easily 
obtained from the relation 

r- = (-i)^exp(y 5 "). ( 4 . 20 ) 

□ 

We would like to mention to conclude this section that it is also possible to 
compute in this framework the matrix elements of more general quasi-local operators, 
i.e of operators obtained as a combination of local spin operators acting on a finite 
number of sites of the lattice. 

5 Correspondence with the results obtained by Alge¬ 
braic Bethe Ansatz 

The XXX chain is a unique example for which the results obtained by SoV can be 
directly compared with the results obtained by ABA. The SU{ 2 ) symmetry of the 
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XXX monodromy matrix provides indeed a correspondence between the chain with 
antiperiodic boundary conditions and the chain with the following twisted boundary 
conditions: 

(5.1) 


0-^+1 = o-ffJi erf, a = x,y,z. 


The eigenstates of this twisted XXX chain can be constructed, in the ABA frame¬ 
work, as the eigenstates of the corresponding twisted transfer matrix. 


r-(A) = tro [ug ro(A)] = A(A) - D{X), 


(5.2) 


by using as usual the operators B{X) (or C'(A)) as creation operators on the reference 
state (j2.14p . More precisely, the state 


R 


i=i 


(5.3) 


is an eigenstate of the transfer matrix 7T(A) if and only if the set of parameters 
{A} = {Ai,..., Xr} satishes the Bethe equations 


R 

n 


Xb- Xg + r] _ a{Xb) 
^b-K-11 d{Xb) ’ 

and the corresponding eigenvalue is 


(5.4) 


r_(A) = a(A) H 11^“^“ + '' 


a=l 


Afe Xg 


a=l 


Afo Aq 


(5.5) 


Let us mention here that there exists an alternative way to construct the eigenstates 
of (j5.2|) . starting instead from the second reference state 




(5.6) 


and using the operators C'(A). Then, if {Ai,..., Xr} satisfies the Bethe equations 
(15.411 . the state 

R 

$({A})) = nC{A,)|0'), (5.7) 

i=i 


is an eigenstate of the transfer matrix corresponding to the eigenvalue 


?_(A) = -r_(A). (5.8) 

Note that the Bethe equations (|5.4I1 issued from the ABA study of the twisted 
chain coincide with the Bethe equations (12.4411 issued from the SoV study of the 
antiperiodic chain, so that the completeness of the former in the ABA framework can 
be derived from the completeness of the latter. It is moreover possible to establish 
an explicit one-to-one correspondence between the eigenstates. 
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To this aim we use the SU{2) symmetry of the XXX monodromy matrix To(A): 
for any U G SU{2), it is easy to see that 


N 

[?7orj^,ro(A)] = 0, where T^ = <8) Un- (5.9) 

n=l 

To establish the relation between two transfer matrices, we consider the following 
unitary matrix: 

= ^ (_1 1 ) • ( 5 . 10 ) 

Using the property Ua^ = a^U it is easy to establish the following similarity trans¬ 
formation between the transfer matrices of two models: 

r-(A) = r^r(A)r-\ (5.11) 

This similarity evidently means that these two transfer matrices are isospectral, and 
that their eigenvectors are in one-to-one correspondence by means of the similarity 
matrix T^. More precisely, we have the following result. 

Theorem 5.1. Let {A} = {Ai,...Ai{} be a solution of the Bethe equations (|5.4I) . 
and let us define the following polynomial of degree R with roots Xj, 1 < j < R: 

R 

Q(A) = n(A-A,). (5.12) 

i=i 

Then the separate state \Q) is an eigenstate of the antiperiodic transfer matrix (|2.5I) 
with eigenvalue 

Q(A) + Q(A) ’ 

whereas the Bethe vector |T({A})) (15.7p is an eigenstate of the twisted transfer 
matrix (I52D with the same eigenvalue (15.1,'ip . These two states are related by means 
of the similarity matrix constructed from ()5.10p as 

\Q) = (-l)^(^-i)2f-^r-^|T({A})). (5.14) 

Proof. For any solution {A} of the Bethe equations we can construct such states 
I Q) and | T({A})). Let us first show that the state | TdA})) is non-trivial. The 
norm of | Q) is given by (I3.62p while for | T({A})) one can use the usual Algebraic 
Bethe Ansatz arguments to obtain a representation in terms of the same Gaudin 
determinant: 


($({A})|$({A})) 



R 


n (Am 

m,n=l 


n (Am 

m^n 


An + V) 
-An) 


det 4 >t-. 
R 


(5.15) 


Since \Q) is a nontrivial state by construction it follows that the corresponding 
Gaudin determinant is non zero and therefore that | T({A})) is non-trivial. 
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The similarity (jS.llI) implies that 


r(A) r-i| T({A})) = r(A) $({A})), (5.16) 


so that r~^| T({A})) is an eigenstate of the antiperiodic transfer matrix T(A) with 
the same eigenvalue t(A) ()5.13|1 as \Q). Hence, since the spectrum of T(A) is 
simple, this state should be proportional to \Q). The proportionality coefficient can 
be computed up to a phase factor from the Gaudin formula. 

It remains to compute this phase factor. Let us consider the simplest separate 
state I 1). It is an eigenstate of T(A) with eigenvalue (i(A)—a(A), and should therefore 
be proportional to r“^| O'). It is easy to see (directly from the construction) that 


( 0 | 1 ) = 1 , 

It leads to a useful relation: 


(o|r-i|o') 



i) = (-A/2)^r;i|o')= S ( 


(5.17) 


(5.18) 


Expressing now | Q) by means of the ABA-type representation (|2.59jl . we obtain 


R R 

rj Q) = (-1)^^ H DiXj) 11) = 2f (-i)M«-i) H (r^ d{x,) t-^) | o' ). 

i=i i=i 

(5.19) 

It is easy to see that 

D{x) r-i = i (a(a) + B{X) + c{x) + d{x)) . (5.20) 

Due to the proportionality between T^ | Q) and | TdA})) only the term containing 
R operators C{Xj) produces a non-zero contribution and we obtain the final result 
for the phase factor. □ 

Hence, Theorem O relates the eigenstates of the antiperiodic transfer matrix, 
which were constructed as separate states in the SoV framework, with the on-shell 
Bethe states for the twisted transfer matrix constructed from ABA (note that Bethe 
equations are crucial here). On the one hand, it provides an easy way to prove 
the completeness of the ABA construction. On the other hand, it gives explicit 
representations valid in the homogeneous limit for the separate states: in particular, 
the equation (|5.18p shows that the state | 1) does not depend on the inhomogeneities. 
A similar result can be obtained for the state | lait), 


|lalt)=SfM ■ (5.21) 

>^=1 V 

This explicit correspondence between these two families of eigenstates can also 
be used to compare the expressions we have just obtained for the form factors of the 
antiperiodic model in the SoV framework with those that can be computed for the 
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twisted model in the ABA framework. Namely, let us fix two given eigenvalues t(A) 
and t'(A) of the antiperiodic (or twisted) transfer matrix, corresponding to two sets 
of solutions {A} = {Ai,..., A/j} and {A'} = {A'^,..., A^,} of the Bethe equations 
(|2.44p with associated Baxter polynomials Qt(A) and Qt'(A). We shall denote as 
usual the corresponding T(A)-eigenstates by | Qr ) andJ^Qr'}, and the corresponding 
71(A)-eigenstates of the form (j5.7p by | 'I'({A})) and | 'I'({A'})). Then, Theorem l5.ll 
implies the following relation between the form factors: 


( Qt I <^n I Qt' 


{®({i)) I < I ®({i'))) - 

($({A})|<|$({A'))) 

2i 

' - { T({A}) I u- I T({A'}) )/2i for R = R' + l, 
< (T({A})|<|T({A'})) ioTR = R', 

_(T({A})|u+|T({A'}))/ 2 i fori?' = ii + l. 


(5.22) 

(5.23) 


Since we are now able to compute these matrix elements independently in the SoV 
and in the ABA framework, and since the obtained representations may appear 
slightly different, it is important to check their effective coincidence. This is explicitly 
done in Appendix [Bj 
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A Proof of Identity [3] 

In this appendix we prove the identity [3] concerning the equality between the quan¬ 
tity and the generalized Slavnov determinant (|3.46l) when {x} is 

a solution to the Bethe equations (I3.29p . We shall proceed in a way similar to the 
proof of the identity [2j 

We consider a set of pairwise distinct Bethe roots {x} = {xi,..., xm}, and a set 
of arbitrary complex numbers {y} = {yi, ... ,yM+s}- As previously we introduce 
the polynomial Q(A) with roots xi,..., xm- 

N 

g(A) = n(A-^i), (A.l) 

i=i 

and the M polynomials Qk{A) obtained from Q(A) as 
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We also introduce polynomials Z(A) of degree S with generic roots zi,..., zs such 
that Za ^ Xj, Za ^ Xj - T], 


s 

Z{\) = X{{\-Za), 

a=l 


ZkW 


Z{\) 

(A - Zk)' 


(A.3) 


As before we consider the auxiliary (2M + S) x (2M + S) matrix C composed of the 
coefficients of the following polynomials; 

2M+S 

Qk{X)Qk{X + v)= E 

i=i 

2M+S 

QkWQ{^ + v)= ^ 

i=i 

2M+S 

Q{X)Q{X + rj)Zk{X)= Y. Ck+ 2 M,jX^-^- (A.4) 

j=i 


Evidently the matrix elements Ckj are zero for k < 2M and j > 2M. However the 
matrix C is invertible and its determinant can easily be computed using the following 
identities: 

2M+S 

^ ^ Cb,a {Xk h) ^b,k Qk{Xk) Qk{Xk 

a=l 

2M-\-S 

^ ^ Cb,a Xf, — ^b,k+M Qixk if) Qk{xk) S^^k Qk{xk) Qk{xk 1]), 

a=l 

2M+S 

^ ^ Cb,a ^k — ^b,k+2M Qi^k ff) Qi^Zk) Zki^Zk) ~\~ ^b,kj 
a=l 

with Ab,k = 0 if 6 > 2M. It means that the product of the matrix C and the 
corresponding Vandermonde matrix is triangular, so that the determinant can be 
easily computed; 


V ({x — r?} U {x} U {z}) det C 

2M+S 

M S 

=n[« {Xk) Qkixk 1]) Q{Xk ??)] Qi^j f]) Zj(^Zj^\ . (A.5) 

k=l j=l 


Hence, 


det C 

2M+S 


, M(M+ 1 ) + S(S- 1 ) 
(- 1 ) 


M 


y{{z}) n [Ql{xk)Qk{xk + rj)] ■ 

k=l 


Let us now compute the following product of two determinants: 


(A.6) 


A 


{x}U{j/} 



det C 

2M+S 


/ g(bi) gik2) \ 

det 2 M +5 I I 

V g(3,i) g(3,2) J 

V({x}U{y}) 


(A.7) 
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As in the case 5 = 0, the M x M block jg given by 
2M+S 

{xk - vT~^) 

a=l 

= Qj{xk) Qjixk + r]) - tJ'E'l^^}{xk) Qj{xk) Qj{xk - r]) 

^jk QkiXk) i^QkiXk 'H) QkiXk • (■^•^) 

Due to the Bethe equations = 0. Another block which turns out to be trivial 

Ji^ 

is the S X M matrix 
2M+S 

G^k^ = Y1 Gj+2M,a (^xl-^ - flE+^{xk) {Xk - vT~^) 
a=l 

= Q{xk) Q{xk + rj) Zj{xk) - Q{xk) Q{xk - v) Zj{xk - f]) 

= 0, (A.9) 

due to the zero factor Q{xk) in both terms. It means that there is no need to 
compute the block as it does not contribute to the determinant. The M x M 

block is diagonal: 

2M+S 

= Y1 Gj+M,a (x^-^ - fiE+^^{xk){Xk - VT~^) 

a=l 

= Qji^k) Qi^k +v) - Q{xk) Qj{xk - v) 

^jk Qfc(^/c) Q{Xk “1“ ^)* (A.10) 

The non-trivial part of this product is contained in the two remaining blocks. The 
(M + S) X M block has a usual form of a Slavnov matrix: 

2M+S 

= Y1 {yt~^ - yE+yivk) ivk - vT~^) 

a=l 

= QjiVk) QjiVk + v)- yE+^^{yk) QjiVk) Qjivk - v) 

=-^Qiyk)Q{yk - ri)'kjk, (A.ii) 

while the (M + 5) x 5 block is 
2M+S 

Gj+2M,a (yt~^ - yE+^ivk) ivk - vT~^) 

a=l 

= Q{yk) {Q{yk + ?/) Zj{yk) - yE^^^{yk) Q{yk - v) Zj{yk - r/)^ . (A.12) 

Now to get rid of the arbitrary polynomials Z(X) we introduce the following S x S 
matrix C defined as 

s 

Zj{X) = Y,Cj,aX"-^- (A.13) 

a=l 
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It is easy to see that 


(A.14) 


detC = (-l)Ty({4). 


Together with the following representation for the block 


QiVk + v) Zjivk) - Qivk - rf) Zj{yk - y) 

S 

= + y) Vk~^ - y^{i}{yk) Qivk - v){yk - yT~^), (A.is) 


a=l 


it leads to the expression (13.471) . 

B Explicit comparison of the form factor representa¬ 
tions issued from SoV and from ABA 

In this appendix we explicitly check that the expressions of the form factors we have 
obtained from our SoV study of the XXX antiperiodic spin chain are consistent with 
the one issued from the ABA study m of the twisted spin chain, i.e. that the 
relation (j5.23p is effectively satisfied. 

The fact that the SoV and the ABA expressions for the form factors coincide if 
\R — R'\ = 1 is quite simple to prove. We shall therefore provide here the details for 
the verification of the last case R = R' only. 

The direct computation of the form factor ( T({A}) | \ 4'({A'})) by ABA leads 

to the following expression: 


(T({A}) I < I T({A'})) = {A'}|{e}) 



(B.l) 


where ^’™^({A}, {A'}|{.^}|^n) is obtained from ^^({A}, {A'}|{^}) (I3.27p by sub¬ 
stituting to A^ in the m-th column of the matrix {A'}|{'C})- The SoV 

computation gives instead: 


(ftk„-|Or.>=A‘’({A).{A')|{{}) 

R 

^ Qrj^'m ~ ^) + dj^'m) Qri^'m + V) M-hm) 

^ n(y -n) ^ 


({A},{A'}|{e}|en). (B.2) 


a(Ad)<3r(Ad - y) 


In the case t(A) = t'(A), the equality between (jB.Ih and (IB.2h is a consequence 
of the following identity: 

O-PAm) Qr(Am T d(Am) Qr(Am T ?/) _ 2a(Am) QriXm T ”^(Am) Qr^Xm') 


®(Am) Qt{X m -y) 


®(Am) Qt{X m -y) 


= 2 . 


(B.3) 


36 















Let us now consider the case R = R' for r(A) / We therefore want to 

compute the difference of the two expressions (jB.ip and (IB.2|1 . 


(Q.|a-|Q,,)-(^({A})|<|$({A'})) 

^ -A a{\'m) [2Qr'(Am -rj)- QT{>^'m - r])] - d{\'^) Qr{\'m + rj) 

a(Am) Qri'^'ra “ V) 

x4''’"AA},{A'}|mi^n), (B.4) 


and show that it vanishes. 

One can first notice that the right hand side of (jB.4p can be seen as the devel¬ 
opment of the determinant of a larger matrix with one more line and column: 


(Q.|a-|g..)-($({A})|<4({A'})) 

^_ detfi+i5(") _ 

F(Ai,...,A^)y(A;,...,Ayna=iQr(AL)’ 


(B.5) 


where 5^”^ is the (i? -|- 1) x {R + 1) matrix of elements, for 1 < j, A: < i?: 


“ QrA^'k) ^Qr,jr(Afc - ??) + j Qr,j{\'k + i 

[S^"^],,R+^=QrMn)QrMn-ri). 

[•?'”’] H+l,i = QriK) {^QA\'k - ») - QriK - ») - +’)>) > 

[S'”L+.,k+. = 


with 7 being an arbitrary complex number. Here we have used the shorthand nota¬ 
tion 

(^-6) 

We now want to rewrite this determinant in terms of the determinant of a {2R + 
1) X {2R + 1) matrix Q = C ■ X. To this aim we introduce the matrices C and X as 
follows. 

The (2i?-|-1) X (2i?-|-1) matrix C is defined by its elements 1 < j, A: < 2R + 1 
such that 


X>~^Cm,j = Qr^mW QT,m.{>^ “ rj), 


i=i 

^ ^ A'^ Cm+R,j — 9m(,^')i 


Mm € {1, 


i=i 

2i?-|-l 

y] A^"^C2ij+i7 


i=i 


h(A), 


Mm E {1,..., i?}, 
.,R], 


(B.7) 

(B.8) 

(B.9) 
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where the functions h{X) and gm{X), 1 < m < R, are given as 

h{X) = {Qr{X) - Qr'{X)) {Qr'{X - v) - Qt{X - ??)) + Qr'{X) Qr'{X - 7 ]), (B.IO) 
2R 

9m{X)= n (B.ll) 

Here gj are general complex numbers which have only to satisfy the condition: 

det [Gj{Xk)]^0, where GkiX) = gkiX) + gk{X + g). (B.12) 

i<j,k<Ra[X) 

In its turn, the (2i? + 1) x {2R + 1) matrix X has for elements 

X^,k = T,iXi), (B.13) 

^j,R+l+k = Tj{Xk) — ^j,2R+l H{Xk), (B.14) 

^j,R+l = + ^j,2R+l{l - (B.15) 

for A: G {1,, i?} and j G {1,..., 2R + 1}, where we have set 

T,{X) = X^-^ + ^{X + gy-\ (B.16) 

HiX) = hiX) + ^hiX + g). (B.17) 

a(A) 

We can write the matrix product of C and X in a block form 

/Gn Gi 2 013 = 0 \ 

C ■ X = Q = j 021 022 023 1 , (B.18) 

\031 032 033 = 0/ 

where we have used that Cm, 2 R+i = 0, for any m G {1,..., 2R}. Here Gii = {Fji^'k)) > 
013 = [Fj{Xk)), 021 = {GjiX'f,)) and G 23 = {Gj{Xk)) are R x R matrices; G 22 = 
{GjiCn)) and 0 i 2 = {Fj{ln)) are Rxl columns; 03 i = (i0(A'^)) and 033 = {K{X'f.)) 
are 1 x R rows; finally 032 = 7 . Here we have dehned 

F,{X) = Q.,,(A) (^Q.,,(A -9) + ^ + 7)) , 

KiX) = QAX) (^Qr'iX - 7 ) - Qr{X -V) + ^ [Qr'A + V) - Qr{X + 7 ?)] ) • 
Note that we have 

H{X'J=K{X'J VmG i?}. (B.19) 

If we now use the block determinant formula to compute the determinant of the 
matrix Q we obtain 

_ det 2 fl+i 0 _ 

V{X,,...,XR)V{X[,...,Xn)ntiQriK) 

= {{Qr\ay\Qr')-{ $ ({A}) I < I $({A'}))) ^^det^^ [G^Xk]], (B. 20 ) 
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once we observe that, thanks to the Bethe equations 


Qr'i^^rn V) ~ Qr'i^rn V): 


one has 

K{X'J = QriX'J (^2QA>^’„^ - r?) - Qr(\'^ - r?) - QAX'^ + r?)) . (B.21) 

This leads to a representation of the quantity ()B.4ll that we want to compute in 
terms of the product of the two determinants of C and <T: 


(Q.|u-|g,.)-($({A})|<|T({A'})) 

__ det2R+i C det2R+i T* _ 

“ B(Ai, ...,Xr) V{X[, ...,X'r) na=i QrA'a) deti<,- [GjiXA] ' 


(B.22) 


We shall now show that 

det T” = 0, (B.23) 

2 _R +1 

which will conclude our proof. In order to do so, let us introduce the (2i?+l) x {2R+ 
1) matrix B with elements given by the coefficients of the following polynomials: 


2R+1 

^ A^"^ BmJ = gr,m(A) Qr.mCA - T]), 

i=i 

2R+1 

Vm E {1,.. 

.,R}, 

(B.24) 

^ A^"^ Bm+R,j = gr,m(A) Qr{X - rj), 
i=i 

2R+1 

Vm E {1,. 


(B.25) 

^ A^-1 B2RJ = Qr{X) Qr{X - V) + HA, 
i=i 

2iJ+l 



(B.26) 

^ X^-^B2R+lJ=h{A- 



(B.27) 

i=i 





Note that this matrix has a non zero determinant: 

/ J?—1 

det B = h{XR) y({Ai}) B({A, + rj}) TT Qr,A^m - r?) ^ 0. (B.28) 

2R+1 ri 

m=l 

Now we can compute the matrix product 

~ /^ll ^12 ^13^“ 

Q = B ■ X = j Q 21 Q 22 Q 2 Z I ) (B.29) 

\^3i ^32 ^33 = 0/ 

where we need to precise only the R x R matrices ^ 23 : 

^^ 23 !, ~ (1 AjR) Ax QT,j(.^j v')i ^ ^ • • •).R}' 

\ / j,k 
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Indeed, the Q blocks are defined above, whereas the Rx R matrix Q 21 and the ii x 1 
column Q 22 have no influence on the computation of the determinant of It is 
then simple to show that, after some row and column exchange, we can write a new 
(2i? + 1) X (2ii + 1) matrix such that 

det G = det G = det A det G23, with G = { .„ ^ ) (B.30) 

2R+1 2R+1 R+2 R-l 023 / 


where ^ is a (ii + 2) x (i? + 2) matrix, A' is a (ii + 2) x (i? — 1) matrix. A” is a 
(ii — 1) X (ii + 2) matrix, and 023 is a (ii — 1) x (ii — 1) diagonal invertible matrix. 
In particular, we have defined 


[023]^ ;, — ^j,kQT,j{^j)QT,j{^j 'n)i Vj)/c G {1, . . . , ii 1}, 
/ 011 012 Orxi\ 


A = 


031 

021 


R 


022 


R 


0 

0 


/ 


(B.31) 

(B.32) 


where ^ 02 i^^ and ^ 022 ^^ are respectively the last row and the last element of 02 i 
and 022- This proves our statement as deti?+ 2 -T = 0. 


References 

[1] H. Bethe. Ziir Theorie der Metalle I. Eigenwerte und Eigenfunktionen Atomkete. 
Zeitschrift fiir Physik, 71:205-226, 1931. 

[2] B. M. McCoy and T. T. Wu. The Two-Dimensional Ising Model. Harvard Univ. 
Press, Cambridge, Massachusetts, 1973. 

[3] R. J. Baxter. Exaetly solved models in statistieal meehanics. Academic Press, 
London, 1982. 

[4] M. Gaudin. La fonction d’onde de Bethe. Masson, 1983. 

[5] L. D. Faddeev and L. A. Takhtajan. Hamiltonian methods in the theory of 
solitons. Springer-Verlag, 1987. 

[6] D. C. Mattis. The Many Body problem: An Encyelopedia of Exactly Solved 
Models in One Dimension. World Scientific, Singapore, 1993. 

[7] B. Sutherland. Beautiful models. World Scientific, 2004. 

[8] L. D. Faddeev. Integrable models in (1 -|- l)-dimensional quantum field theory. 
In J. B. Zuber and R. Stora, editors, Les Houehes 1982, Recent advanees in field 
theory and statistieal meehanics, pages 561-608. Elsevier Science Publ., 1984. 

[9] L. D. Faddeev. How algebraic Bethe ansatz works for integrable model, hep- 
th/9605187, 1996. 

[10] Beisert N. et al. Review of AdS/CFT Integrability; An Overview. Lett. Math. 
Phys., 3:99, 2012. 


40 



[11] E. K. Sklyanin and L. D. Faddeev. Quantum mechanical approach to completely 
integrable field theory models. Sov. Phys. Dokl.^ 23:902-904, 1978. 

[12] L. D. Faddeev, E. K. Sklyanin, and L. A. Takhtajan. Quantum inverse problem 
method I. Theor. Math. Phys., 40:688-706, 1979. Translated from Teor. Mat. 
Fiz. 40 (1979) 194-220. 

[13] E. K. Sklyanin. Method of the inverse scattering problem and the non-linear 
quantum Schrodinger equation. Sov. Phys. Dokl, 24:107-109, 1979. 

[14] P. P. Kulish and E. K. Sklyanin. Quantum spectral transform method. Recent 
developments. Lectures Notes in Physics, 151:61-119, 1982. 

[15] C. N. Yang. Some exact results for the many-body problem in one dimension 
with repulsive delta-function interaction. Phys. Rev. Lett., 19:1312-1314, 1967. 

[16] R. J. Baxter. One-dimensional anisotropic Heisenberg chain. Phys. Rev. Lett., 
26:834, 1971. 

[17] R. J. Baxter. Eight-vertex model in lattice statistics. Phys. Rev. Lett., 26:832- 
833, 1971. 

[18] M. Jimbo. A q-difference analogue of U{g) and the Yang-Baxter equation. Lett. 
Math. Phys., 10:63-69, 1985. 

[19] M. Jimbo. A q-analogue of U{gl{N -\- 1)), Hecke algebra, and the Yang-Baxter 
equation. Lett. Math. Phys., 11:247-252, 1986. 

[20] P. P. Kulish and N. Yu. Reshetikhin. Quantum linear problem for the sine- 
Gordon equation and higher representations. Zap. Nauch. Sem. LOMI, 101:101- 
110, 1981. Translation in J. Sov. Math. 23 (1983) 2435-41. 

[21] V. G. Drinfel’d. Hopf algebras and the quantum Yang-Baxter equation. Soviet 
Math. Dokl, 32(l):254-258, 1985. 

[22] V. G. Drinfel’d. Quantum groups. In Proc. Internat. Congress of Math., Berke¬ 
ley, USA, 1986, pages 798-820. AMS, 1987. 

[23] E. K. Sklyanin. The quantum Toda chain. Lectures Notes in Physics, 226:196- 
233, 1985. 

[24] E. K. Sklyanin. Separation of variables. New trends. Prog. Theor. Phys., 118:35- 
60, 1995. arXiv:solv-int/9504001. 

[25] E. K. Sklyanin. Functional Bethe Ansatz. In B.A. Kupershmidt, editor, In¬ 
tegrable and Superintegrable Systems, pages 8-33. World Scientific, Singapore, 
1990. 

[26] E. K. Sklyanin. Quantum inverse scattering method. Selected topics. In Mo- 
Lin Ge, editor. Quantum Group and Quantum Integrable Systems, pages 63-97. 
Nankai Lectures in Mathematical Physics, World Scientific, 1992. arXiv:hep- 
th/9211111. 


41 



[27] S. Niekamp, T. Wirth, and H. Frahm. The XXZ model with anti-periodic 
twisted boundary conditions. J. Phys. A: Math. Theor., 42:195008, 2009. 

[28] H. Frahm, A. Seel, and T. Wirth. Separation of variables in the open XXX 
chain. Nucl. Phys. B, 802:351, 2008. 

[29] O. Babelon, D. Bernard, and F. A. Smirnov. Quantization of solitons and the 
restricted sine-Gordon model. Commun.Math.Phys., 182:319-354, 1996. 

[30] O. Babelon, D. Bernard, and F. A. Smirnov. Null-vectors in integrable field 
theory. Commun.Math.Phys., 186:601, 1997. 

[31] F. A. Smirnov. Quasi-classical study of form factors in hnite volume, hep- 
th/9802132, 1998. 

[32] O. Babelon. Universal exchange algebra for Bloch waves and Liouville theory. 
J. Phys. A, 37:303, 2004. 

[33] G. von Gehlen, N. lorgov, S. Pakuliak, and V. Shadura. The Baxter-Bazhanov- 
Stroganov model: separation of variables and the Baxter equation. J. Phys. A: 
Math. Theor., 39:7257, 2006. 

[34] G. von Gehlen, N. lorgov, S. Pakuliak, V. Shadura, and Yu. Tykhyy. Form- 
factors in the Baxter-Bazhanov-Stroganov model I: norms and matrix elements. 
J. Phys. A: Math. Theor., 40:14117, 2007. 

[35] G. von Gehlen, N. lorgov, S. Pakuliak, V. Shadura, and Yu. Tykhyy. Form- 
factors in the Baxter-Bazhanov-Stroganov model II: Ising model on the finite 
lattice. ,J. Phys. A: Math. Theor., 41:095003, 2008. 

[36] G. von Gehlen, N. lorgov, S. Pakuliak, and V. Shadura. Factorized finite-size 
Ising model spin matrix elements from separation of variables. J. Phys. A: 
Math. Theor., 42:304026, 2009. 

[37] S. E. Derkachov, G. P. Korchemsky, and A. N. Manashov. Baxter Q-operator 
and separation of variables for the open SL(2,M) spin chain. JHEP, 10:053, 
2003. 

[38] S. E. Derkachov, G. P. Korchemsky, , and A. N. Manashov. Separation of 
variables for the quantum SL(2,M) spin chain. JHEP, 07:047, 2003. 

[39] S. E. Derkachov, G.P. Korchemsky, and A. N. Manashov. Noncompact Heisen¬ 
berg spin magnets from high-energy QGD. I. Baxter Q-operator and separation 
of variables. Nucl. Phys. B, 617:375-440, 2001. 

[40] A. Bytsko and J. Teschner. Quantization of models with non-compact quantum 
group symmetry. Modular XXZ magnet and lattice sinh-Gordon model. J. 
Phys. A, 39:12927, 2006. 

[41] G. Niccoli and J. Teschner. The Sine-Gordon model revisited I. J. Stat. Mech., 
page P09014, 2010. 


42 



[42] G. Niccoli. Reconstruction of Baxter Q-operator from Sklyanin SOV for cyclic 
representations of integrable quantum models. Nucl. Phys. B, 835:263-283, 
2010. 

[43] G. Niccoli. Completeness of Bethe Ansatz by sklyanin SOV for cyclic represen¬ 
tations of integrable quantum models. JHEP, 03:123, 2011. 

[44] G. Niccoli. Non-diagonal open spin-1/2 XXZ quantum chains by separation of 
variables: Complete spectrum and matrix elements of some quasi-local opera¬ 
tors. J. Stat. Mech., page P10025, 2012. 

[45] N. Grosjean and G. Niccoli. The r 2 -model and the chiral Potts model revisited: 
completeness of Bethe equations from Sklyanin’s SOV method. J. Stat. Mech., 
page P11005, 2012. 

[46] G. Niccoli. Antiperiodic spin-1/2 XXZ quantum chains by separation of vari¬ 
ables: Complete spectrum and form factors. Nucl. Phys. B, 870:397-420, 2013. 
arXiv: 1205.4537. 

[47] G. Niccoli. An antiperiodic dynamical six-vertex model: I. Complete spectrum 
by SOV, matrix elements of the identity on separate states and connections to 
the periodic eight-vertex model. J. Phys. A: Math. Theor., 46:075003, 2013. 
arXiv:1207.1928. 

[48] G. Niccoli. Form factors and complete spectrum of XXX antiperiodic higher 
spin chains by quantum separation of variables. J. Math. Phys., page 053516, 
2013. 

[49] S. Faldella and G. Niccoli. SOV approach for integrable quantum models asso¬ 
ciated with general representations on spin-1/2 chains of the 8-vertex reflection 
algebra. J. Phys. A: Math. Theor., 47:115202, 2014. 

[50] S. Faldella, N. Kitanine, and G. Niccoli. Complete spectrum and scalar products 
for the open spin-1/2 XXZ quantum chains with non-diagonal boundary terms. 
J. Stat. Mech., page POlOll, 2014. 

[51] G. Niccoli and V. Terras. Antiperiodic XXZ chains with arbitrary spins: 
Complete eigenstate construction by functional equations in separation of 
variables. Lett. Math. Phys., pages DOI 10.1007/sl 1005-015-0759-9, 2015. 
arXiv:1411.6488. 

[52] V. Tarasov and A. Varchenko. Completeness of Bethe vectors and differ¬ 
ence equations with regular singular points. Int. Math. Research Notices, 
1995(13) :637-669, 1995. 

[53] E. Mukhin, V. O. Tarasov, and A. Varchenko. Bethe algebra of homogeneous 
XXX Heisenberg model has simple spectrum. Comm. Math. Phys., 288(l):l-42, 
2009. 

[54] F. A. Smirnov. Form factors in completely integrable models of quantum field 
theory. World Scientific, Singapore, 1992. 


43 



[55] V. E. Korepin, N. M. Bogoliubov, and A. G. Izergin. Quantum inverse scattering 
method and correlation functions. Cambridge University Press, 1993. 

[56] M. Jimbo and T. Miwa. Algebraic analysis of solvable lattice models. Number 85 
in CBMS Regional Conference Series in Mathematics. AMS, Providence, RI, 
1995. 

[57] N. Kitanine, J. M. Maillet, and V. Terras. Form factors of the XXZ Heisenberg 
spin-1/2 finite chain. Nucl. Phys. B, 554:647-678, 1999. 

[58] N. Kitanine, J. M. Maillet, and V. Terras. Correlation functions of the XXZ 
Heisenberg spin-1/2 chain in a magnetic field. Nucl. Phys. B, 567:554-582, 2000. 

[59] N. Kitanine, J. M. Maillet, N. A. Slavnov, and V. Terras. Spin-spin correlation 
functions of the XXZ-1/2 Heisenberg chain in a magnetic field. Nucl. Phys. B, 
641:487-518, 2002. 

[60] N. Kitanine, J. M. Maillet, N. A. Slavnov, and V. Terras. Dynamical correlation 
functions of the XXZ spin-1/2 chain. Nucl. Phys. B, 729:558-580, 2005. hep- 
th/0407108. 

[61] H. Boos, M. Jimbo, T. Miwa, F. Smirnov, and Y. Takeyama. Hidden Grassmann 
structure in the XXZ model. Comm. Math. Phys., 272:263-281, 2007. 

[62] N. Kitanine, K. K. Kozlowski, J. M. Maillet, G. Niccoli, N. A. Slavnov, and 
V. Terras. Correlation functions of the open XXZ chain: I. J. Stat. Mech. 
Theory Exp., page P10009, 2007. 

[63] N. Kitanine, K. K. Kozlowski, J. M. Maillet, G. Niccoli, N. A. Slavnov, and 
V. Terras. Correlation functions of the open XXZ chain: H. J. Stat. Mech. 
Theory Exp., page P07010, 2008. 

[64] H. Boos, M. Jimbo, T. Miwa, F. Smirnov, and Y. Takeyama. Hidden Grassmann 
Structure in the XXZ Model H: Creation Operators. Comm. Math. Phys., 
286:875-932, 2009. 

[65] J. M. Maillet and V. Terras. On the quantum inverse scattering problem. Nucl. 
Phys. B, 575:627-644, 2000. hep-th/9911030. 

[66] A. G. Izergin, N. Kitanine, J. M. Maillet, and V. Terras. Spontaneous mag¬ 
netization of the XXZ Heisenberg spin-1/2 chain. Nucl. Phys. B, 554:679-696, 
1999. solv-int/9812021. 

[67] N. Kitanine, J. M. Maillet, N. A. Slavnov, and V. Terras. Master equation for 
spin-spin correlation functions of the XXZ chain. Nucl. Phys. B, 712:600-622, 
2005. hep-th/0406190. 

[68] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, and V. Terras. Alge¬ 
braic Bethe ansatz approach to the asymptotic behavior of correlation functions. 
J. Stat. Mech. Theory Exp., page P04003, 2009. 


44 



[69] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, and V. Terras. On 
the thermodynamic limit of form factors in the massless XXZ Heisenberg chain. 
J. Math. Phys., 50:095209, 2009. 

[70] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, and V. Terras. A 
form factor approach to the asymptotic behavior of correlation functions in crit¬ 
ical models. J. Stat. Mech. Theory Exp., page P12010, 2011. arXiv:1110.0803. 

[71] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, and V. Terras. 
The thermodynamic limit of particle-hole form factors in the massless XXZ 
Heisenberg chain. J. Stat. Mech.: Theory Exp., page P05028, 2011. 

[72] N. Kitanine, K. K. Kozlowski, J. M. Maillet, and V. Terras. Large-distance 
asymptotic behaviour of multi-point correlation functions in massless quantum 
models. J. Stat. Mech., page P05011, 2014. 

[73] K. K. Kozlowski and J. M. Maillet. Microscopic approach to a class of Id 
quantum critical models. arXiv;1501.07711, 2015. 

[74] K. K. Kozlowski, J. M. Maillet, and N. A. Slavnov. Long-distance behavior of 
temperature correlation functions of the quantum one-dimensional Bose gas. J. 
Stat. Mech.: Theory Exp., page P03018, 2011. arXiv;1011.3149. 

[75] K. K. Kozlowski, J. M. Maillet, and N. A. Slavnov. Correlation functions for 
one-dimensional bosons at low temperature. J. Stat. Mech.: Theory Exp., page 
P03019, 2011. 

[76] K. K. Kozlowski and V. Terras. Long-time and large-distance asymptotic be¬ 
havior of the current-current correlators in the non-linear Schrodinger model. 
J. Stat. Mech.: Theory Exp., page P09013, 2011. arXiv; 1101.0844. 

[77] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, and V. Terras. 
Form factor approach to dynamical correlation functions in critical models. J. 
Stat. Mech. Theory Exp., page P09001, 2012. arXiv:1206.2630. 

[78] M. Dugave, F. Gohmann, and K. K. Kozlowski. Thermal form factors of the 
XXZ chain and the large-distance asymptotics of its temperature dependent 
correlation functions. J. Stat. Mech., page P07010, 2013. 

[79] M. Dugave, F. Gohmann, and K. K. Kozlowski. Low-temperature large-distance 
asymptotics of the transversal two-point functions of the XXZ chain. J. Stat. 
Mech., page P04012, 2014. 

[80] J. S. Caux, R. Hagemans, and J. M. Maillet. Gomputation of dynamical corre¬ 
lation functions of Heisenberg chains: the gapless anisotropic regime. J. Stat. 
Mech. Theory Exp., page P09003, 2005. 

[81] J. S. Gaux and J. M. Maillet. Computation of dynamical correlation functions 
of Heisenberg chains in a magnetic field. Phys. Rev. Lett., 95:077201, 2005. 

[82] A. G. Izergin. Partition function of the six-vertex model in a finite volume. Sov. 
Phys. Dokl., 32:878-879, 1987. 


45 



[83] N. A. Slavnov. Calculation of scalar products of wave functions and form factors 
in the framework of the algebraic Bethe Ansatz. Theor. Math. Phys., 79:502- 
508, 1989. 

[84] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, and V. Terras. 
Riemann-Hilbert approach to a generalized sine kernel and applications. Comm. 
Math. Phys., 291:691-761, 2009. 

[85] N. Kitanine, K. K. Kozlowski, J. M. Maillet, N. A. Slavnov, and V. Terras. 
Thermodynamic limit of particle-hole form factors in the massless XXZ Heisen¬ 
berg chain. arXiv: 1003.4557. 

[86] N. Grosjean, J. M. Maillet, and G. Niccoli. On the form factors of local operators 
in the lattice sine-Gordon model. J. Stat. Mech.: Theory Exp., page P10006, 
2012. 

[87] N. Grosjean, J. M. Maillet, and G. Niccoli. On the form factors of local operators 
in the Bazhanov-Stroganov and chiral Potts models. Annales Henri Poincare, 
16:1103, 2015. 

[88] D. Levy-Bencheton, G. Niccoli, and V. Terras. Antiperiodic dynamical 6-vertex 
model by separation of variables II: Functional equation and form factors, in 
preparation. 

[89] I. Kostov. Glassical limit of the three-point function of n= 4 supersymmetric 
yang-mills theory from integrability. Phys. Rev. Letters, 26:261604, 2012. 

[90] I. Kostov. Three-point function of semiclassical states at weak coupling. J. 
Phys. A: Math. Theor., 45:494018, 2012. 

[91] I. Kostov and Y. Matsuo. Inner products of Bethe states as partial domain wall 
partition function. JHEP, 10:168, 2012. 

[92] O. Foda and M. Wheeler. Variations on Slavnov’s scalar product. JHEP, 10:096, 

2012. 

[93] A. Garbali. The scalar product of XXZ spin chain revisited. Application to the 
ground state at A = 1/2. arXiv:1411.2938. 

[94] V. Terras. Drinfel’d twists and functional Bethe Ansatz. Lett. Math. Phys., 
48:263-276, 1999. math-ph/9902009. 

[95] J. M. Maillet and J. Sanchez de Santos. Drinfel’d twists and algebraic Bethe 
Ansatz. In L. D. Faddeev’s Seminar on Mathematical Physics, pages 137-178. 
Amer. Math. Soc. Transl. Ser. 2, 201, Amer. Math. Soc., Providence, RI, 2000. 
q-alg/9612012. 

[96] N. Kitanine, J. M. Maillet, and G. Niccoli. Open spin chains with generic 
integrable boundaries: Baxter equation and Bethe ansatz completeness from 
separation of variables. J. Stat. Mech., page P05015, 2014. 

[97] E. Bettelheim and I. Kostov. Semi-classical analysis of the inner product of 
Bethe states. J. Phys A: Math. Theor., 24:245401, 2014. 


46 



